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A  deterministic  model  of  succession  in  aspen  forests 
was  developed  using  existing  data  and  intuition.  The 
degree  of  uncertainty,  which  was  determined  by  allow- 
ing the  parameter  values  to  vary  at  random  within 
limits,  was  larger  than  desired.  Analysis  of  model  sen- 
sitivity to  changes  in  parameter  values  was  made  and 
the  results  presented.  These  results  have  indicated 
areas  of  needed  research.  The  model  responds 
realistically  to  various  management  techniques  and 
could  be  an  aid  to  resource  managers  in  their  decision 
making  process. 
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Aspen  Succession  in  the 
Intermountain  West:  A 
Deterministic  Model 

Dale  L  Bartos,  Frederick  R.  Ward, 
and  George  S.  Innis 


INTRODUCTION 

The  aspen  (Populus  tremuloides  Michx.)  ecosystem  is 
prevalent  in  the  mountainous  West  (Little  1971).  This  system 
produces  multiple  resources  including:  water,  wildlife  habitat, 
recreational  sites,  summer  grazing  for  domestic  livestock,  wood 
fiber,  and  esthetics. 

To  manage  the  aspen  system  and  its  resources  intelligently, 
we  must  understand  the  natural  forces  acting  upon  and  within 
the  ecosystem.  When  we  understand  the  dynamics  of  the  vege- 
tation, we  can  make  valid  judgments  concerning  management 
of  these  lands  for  the  resources  mentioned  above.  The  process 
of  succession  also  needs  to  be  better  understood  so  that  the 
successional  position  of  a  particular  stand  can  be  identified  and 
the  stand  managed  accordingly. 

Much  of  successional  theory  has  been  developed  for  mesic, 
temperate  forests.  Many  of  these  theories  have  been  incor- 
porated into  simulation  models  tracking  the  dynamics  of  the 
forested  system  following  disturbance.  Shugart  and  West  (1980) 
review  much  of  this  literature  and  Trimble  and  Shriner  (1981) 
inventory  many  of  the  published  and  unpublished  models. 

More  than  200  years  is  required  for  conifers  to  invade  and 
dominate  aspen  forests.  Long-term  intensive  studies  are  not 
feasible  to  answer  pressing  and  current  management  problems. 
Using  existing  data  and  intuition,  Bartos  (1973,  1978)  developed 
a  simulation  model  as  an  aid  to  understanding  the  dynamics  of 
the  aspen  system.  We  have  revised  this  model  by  reducing  the 
number  of  parameters,  while  still  retaining  the  structure  and 
behavior  of  the  model. 

The  model  presented  here  considers  a  site  as  a  homogeneous 
unit  on  which  plants  (trees,  shrubs,  and  herbs)  are  growing. 
Because  the  model  addresses  average  conditions  on  an  area,  it 
is  more  like  the  model  of  Noble  and  Slatyer  (1977)  than  the 
family  of  models  produced  by  Botkin,  his  associates  and 
followers  (Botkin  1977;  Botkin  and  others  1972a,  b;  Shugart 
and  West  1977).  Our  model  is  concerned  with  the  alternative 
fates  of  an  aspen/conifer  site  following  a  major  disturbance 
(fire  or  clearcut).  These  alternatives,  described  by  Mueggler 
(1976),  include  succession  to  conifers,  remaining  in  aspen,  or 
succession  to  meadow,  again  suggesting  the  Noble  and  Slatyer 
model.  To  date,  our  model  treats  only  the  first  alternative.  Im- 
plementation of  other  alternatives  is  straightforward,  depending 
only  on  the  identification  of  mechanisms  responsible  for  the 
alternate  pathways. 

Our  model  is  most  like  that  of  Sperger  (1980),  although  the 
latter  model  is  very  elaborate.  Sperger  includes  many  more  size 


classes  in  his  model  so  that  he  can  predict  the  size  class  distri- 
bution of  the  forest  at  any  age.  We  deal  with  aspen  and  conifer 
reproduction  in  three  different  size  classes,  and  beyond  that  we 
consider  only  biomass.  Therefore,  our  model  requires  less  data 
than  Sperger's  model  for  a  particular  run. 

With  all  the  published  successional  models,  why  do  we  pre- 
sume to  build  another?  First,  much  of  the  utility  in  a  modeling 
exercise  derives  from  the  design  and  construction  of  the  model 
(Innis  1972;  Odum  1981).  This  utility  is  denied  one  who  at- 
tempts to  use  another's  model.  Because  models  are  specific  to 
their  applications,  model  building  begins  with  statements  of  the 
objectives.  Each  of  the  published  models  has  different  objec- 
tives and,  consequently,  different  structure.  We  must  appreciate 
this  dependence  on  objectives  when  evaluating  a  model. 

In  addition  to  the  "standard"  run  showing  succession  in  the 
aspen  forest  following  a  disturbance,  this  paper  includes 
responses  of  the  model  to  common  management  techniques.  As 
an  aid  to  interpretation  of  model  results,  an  indication  of  the 
distribution  of  the  model  results  as  a  function  of  uncertainty  of 
the  parameter  values  used  is  presented.  The  model  is  very  sensi- 
tive to  changes  in  some  parameter  values,  and  this  sensitivity  is 
discussed. 

Model  Objective 

Our  objective  was  to  develop  a  streamlined  simulation  model 
that  would  predict  the  dynamic  nature  of  vegetation  compo- 
nents within  the  aspen  system  and  how  these  components 
change  during  the  successional  process.  This  model  should  ap- 
ply to  most  aspen  lands  in  the  western  United  States  and 
should  aid  natural  resource  managers  in  their  decisionmaking 
process.  In  addition,  a  good  responsive  model  would  aid  in 
studying  the  aspen  system.  Such  a  model  would  provide  an  ex- 
cellent basis  on  which  to  develop  a  problem  analysis  for  future 
research.  Individual  studies  could  be  conducted  to  more  accu- 
rately estimate  the  most  sensitive  parameters. 

Several  terms  in  this  statement  of  objectives  require  further 
elaboration:  1.  "Dynamic"  implies  changes  over  time  in  the 
various  vegetation  components  of  the  system.  2.  "Vegetation 
component"  refers  to  only  the  aboveground  biomass  for  wood 
(aspen  and  conifers)  and  herbage  (shrub  and  herb  species).  3. 
"Successional  process"  is  the  dynamics  of  a  vegetation  unit  as 
it  progresses  towards  a  community  that  is  in  a  steady  state  (cli- 
max) relationship  with  the  environment.  We  are  concerned  only 
with  secondary  succession — where  the  soil  may  remain  at  a  fair- 
ly advanced  stage  of  development  while  the  plant  community  is 
set  back  to  an  earlier  stage  because  of  a  disturbance. 

EXCHANGE  Rec'd 


4.  "General  model"  implies  flexibility  to  apply  the  model  to 
different  sites  which  have  similar  environmental  conditions. 
The  model  needs  to  assimilate  disconnected  bits  of  information 
and  predict  within  reason  the  results  of  management  practices. 

While  this  objective  statement  sets  the  general  tone  of  the 
modeling  effort,  a  more  explicit  set  of  issues  must  be 
addressed: 

1.  Does  the  modeled  successional  pattern  fit  the  intuition  of 
experts? 

2.  Does  the  model  respond  reasonably  to  management 
actions? 

3.  Without  major  influences  from  conifers  and  herbs,  is  an 
aspen  community  able  to  stabilize  for  long  periods  of  time? 

4.  What  are  the  sensitivities  of  these  predictions  to 
parameter  estimates,  site  condition  assignments,  initial  condi- 
tions, and  functional  forms  of  flows? 

5.  What  additional  research  is  needed  to  make  this  model  a 
more  accurate  simulation  of  the  consequences  of  succession? 

Assumptions 

Several  assumptions  were  made  to  facilitate  initial  model 
development  and  subsequent  updates.  The  initial  assumptions 
are:  1.  The  system  stressed  will  have  conifers  on  or  near  it  to 
provide  a  seed  source  when  the  site  is  disturbed.  2.  Aspen  are 
present;  therefore,  a  root  source  exists  to  provide  a  flush  of 
aspen  suckers  when  the  site  is  disturbed.  3.  The  damaged  site 
will  be  small  enough  to  be  homogeneous  with  respect  to  soil 


characteristics  and  abiotic  factors.  4.  Disturbance  will  be  by 
burning;  some  root  crowns  of  perennial  herbs  will  not  be 
killed  and  total  herb  production  will  be  temporarily  stimu  - 
lated  when  the  site  is  disturbed.  (The  community  could  quite 
conceivably  be  stressed  by  cutting  the  overstory  trees  and 
similar  responses  might  be  observed.) 

MODEL  STRUCTURE 

The  major  components  of  interest  in  this  model  are  aspen, 
conifers,  shrubs,  and  herbs.  These  regenerate,  grow,  and  die  at 
rates  affected  by  the  various  components  (fig.  1).  Soil  and 
climate  are  not  explicitly  treated  in  the  model  but  are  reflected 
in  site-specific  parameters  such  as  growth  rate,  regeneration 
rate,  and  maximum  biomass  of  the  various  compartments.  The 
mathematical  model  is  a  set  of  simultaneous  difference  equa- 
tions. The  model  was  coded  in  FORTRAN  V  and  runs  on  a 
variety  of  computers.  (A  complete  listing  of  the  computer  program 
is  presented  in  appendix  A.)  The  model  is  deterministic,  not 
stochastic,  in  nature.  Parameter  values  used  (see  appendix  C) 
reflect  long-term  averages  adjusted  for  the  step  size,  which  we 
set  at  1  year. 

Primary  concern  is  aboveground  biomass  (kg/ha)  of  the  four 
major  state  variables.  Aspen  suckers  and  conifer  seedlings, 
however,  are  treated  in  three  size  categories  on  the  basis  of 
density  because  growth,  mortality,  and  interactions  with  other 
compartments  are  very  different  for  the  suckers  and  seedlings 
than  for  mature  individuals  and  should  be  handled  in  a  differ- 
ent fashion. 


Source  or  Sink 


Figure  1.— Summary  flow  diagram  of  major  compartments  in  the  ASPEN  model. 
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The  aspen  reproduction  categories  are:  ASP1,  aspen  suckers 
less  than  0.5  meters  tall;  ASP2,  suckers  between  0.5  and  2 
meters  tall;  and  ASP3,  suckers  over  2  meters  tall  but  less  than 
5  cm  d.b.h.  All  aspen  reproduction  is  considered  to  be 
vegetative.  For  conifers  the  three  categories  of  seedlings  are 
CON1,  CON2,  and  CON3,  defined  according  to  the  same  size 
structure  as  the  aspen  suckers.  CON1  seedlings  are  considered 
to  be  established  (to  have  persisted  on  the  site  for  at  least  3 
years)  to  avoid  dividing  them  into  several  subgroups.  Time  for 
passage  through  these  subgroups  is  accounted  for  by  a  time 
delay.  Thus,  no  regeneration  into  CON1  is  allowed  for  a  few 
years  (3  to  10  depending  on  the  site)  after  the  initial  dis- 
turbance. The  numbers  of  individuals  graduating  from  ASP3 
and  CON3  are  converted  to  biomass  using  average  weights  of 
5  cm  d.b.h.  aspen  and  conifer  and  constitute  the  regenera- 
tion flows  from  the  source  to  the  biomass  compartments  of 
aspen  and  conifer.  These  two  major  compartments  contain  the 
biomass  only  of  individuals  larger  than  5  cm  d.b.h. 

In  the  case  of  shrubs,  mature  plants  were  not  considered  suf- 
ficiently different  from  the  young  ones  to  warrant  separate 
treatments;  so  the  shrub  compartment  contains  individuals  of 
all  sizes. 

The  herb  compartment  was  subdivided  into  annuals  and  per- 
ennials. Because  only  the  aboveground  biomass  is  considered, 
these  compartments  are  totally  depleted  each  year.  Production 
was  not  separated  into  reproduction  and  growth. 

Flows  representing  regeneration,  growth,  production,  or 
graduation  are  controlled  by  the  product  of  a  maximum  rate 
with  two  or  more  inhibiting  functions  having  values  which 
range  over  the  interval  from  0  to  1 .  For  example,  the  gradua- 
tion of  conifer  seedlings  from  CON3  to  CON  is  given  by 
FLOW27  =  CON3  *  P(101)  *  Zl  *  Z2  where  P(101)  gives  the 
maximum  graduation  rate  allowed,  Zl  =  1  -  0.5  *  PASP  gives 
the  restriction  of  the  aspen,  and  Z2  =  1  -  0.7  *  PCON  gives 
the  restriction  of  the  conifers.  The  maximum  rates  are  given  by 
parameters  reflecting  site  characteristics,  and  the  functions 
represent  restrictions  imposed  on  that  flow  by  various  compart- 
ments. The  independent  variable  for  each  restricting  function  is 
taken  to  be  the  ratio  of  the  present  value  to  the  maximum  pos- 
sible value  of  the  inhibiting  compartments,  rendering  these 
functions  independent  of  site. 

With  one  exception,  the  influence  of  aspen  on  conifer  repro- 
duction which  will  be  discussed  later,  these  inhibiting  functions 
are  nonincreasing.  When  the  inhibiting  compartment  has  a 
value  of  zero,  the  ratio  of  present  to  maximum  value  is  zero 
and  the  inhibiting  function  has  a  value  of  1  (no  restriction).  As 
the  value  of  the  inhibiting  compartment  increases,  the  ratio  of 
that  value  to  the  maximum  increases  and  the  value  of  the  in- 
hibiting function  decreases  below  1  so  that  there  is  a  more 
stringent  restriction. 

The  inhibiting  functions  were  suggested  in  graphic  form  by 
members  of  the  Aspen  Ecosystem  Project  (INT-RWU-1751)  at 
the  Intermountain  Station's  Forestry  Sciences  Laboratory  in 
Logan,  Utah.  We  then  fitted  these  curves  with  functions  and 
they  were  approved  by  consensus  of  the  project  scientists  as  ex- 
hibiting appropriate  behavior. 

All  compartments  other  than  sources  and  sinks  were  consid- 
ered to  have  potential  influence  on  a  flow.  If  this  influence  was 
considered  quite  minor,  as  in  the  case  of  herbage  on  growth  of 
large  conifers,  it  was  omitted.  In  cases  where  intensive  in- 
vestigation has  yielded  more  detailed  knowledge,  the  influences 
were  considered  separately,  but  where  less  was  known  the  in- 
fluences were  pooled.  For  instance,  the  three  aspen  sucker  com- 


partments separately  affect  the  regeneration  of  ASP1,  but  the 
three  conifer  seedling  compartments  are  pooled  into  one  com- 
partment, which  has  a  single  restriction  on  ASP1  regeneration. 
The  detailed  discussion  of  the  model  structure  includes  graphs 
of  the  functions  restricting  the  flows  into  the  various  compart- 
ments. Each  graph  is  accompanied  by  its  equation,  and  the 
flow  equations  into  and  out  of  each  compartment  are  given. 

Each  mortality  flow  in  the  model  is  a  fixed  percentage  of  the 
value  of  the  compartment  from  which  the  flow  originates;  for 
example,  the  mortality  of  conifers  is  given  by  the  equation 
FLOW6  =  CON  *  P(105).  In  the  biomass  compartments 
"mortality"  is  interpreted  to  include  the  loss  of  portions  of  a 
plant  as  well  as  of  entire  plants.  These  mortality  rates  are  site- 
specific  parameters.  During  the  model  building  process  there 
was  much  discussion  concerning  the  mechanisms  which  limit  a 
compartment  to  some  maximum  value.  In  the  absence  of  firm 
knowledge  as  to  whether  (1)  the  mortality  rate  increases  to  sur- 
pass the  growth  rate,  (2)  the  growth  rate  drops  to  the  mortality 
rate,  or  (3)  the  growth  rate  decreases  over  time  as  the  mortality 
rate  increases,  it  was  decided  to  hold  mortality  constant  and 
decrease  growth  from  some  maximum  value  by  means  of  the 
various  restrictions.  In  the  aspen  compartment,  for  example, 
maximum  growth  rate  is  set  at  8  percent  per  annum  and  mor- 
tality rate  is  set  at  3  percent  per  annum.  This  allows  for  a  net 
growth  of  8  percent  -  3  percent  =  5  percent  per  annum  in  the 
extreme  case  that  little  aspen  and  no  conifer  is  present  to  curtail 
the  growth  rate.  A  heavy  conifer  stand,  however,  would  restrict 
aspen  growth  rate  to  nearly  zero,  yielding  a  net  growth  of 
0  percent  -  3  percent  =  -  3  percent  per  annum,  the  mortality 
rate. 

Large-scale  losses,  such  as  those  induced  by  disease  or  insect 
infestation,  are  outside  the  purview  of  the  model.  One  could 
treat  such  high  mortality,  however,  by  stopping  the  simulation 
and  restarting  with  new  initial  values  that  reflect  the  appro- 
priate loss. 

At  the  beginning  of  a  simulation  run  the  computer  assigns 
and  prints  the  maximum  possible  values  of  various  compart- 
ments (code  1 1 100-13700).  The  maximum  value  of  a  compart- 
ment represents  the  level  that  compartment  could  contain  in  the 
absence  of  any  competition  from  other  compartments.  Input 
parameters  give  the  maximum  values  of  mature  aspen,  mature 
conifer,  shrubs,  annuals,  perennials,  the  three  aspen  sucker 
classes,  the  three  conifer  seedling  categories,  the  pooled 
biomass  equivalent  of  the  conifer  seedlings,  the  pooled  biomass 
equivalent  of  the  aspen  suckers,  and  the  combined  biomass 
possible  for  shrubs  and  aspen  suckers.  Maximum  herbage  is  the 
sum  of  the  maxima  of  the  annual  and  perennial  compartments. 

Some  initial  values  are  assigned  through  parameters  and 
others  are  computed  from  those  values  assigned  (code 
14400-18100).  Those  assigned  are  for  aspen  and  its  three  sucker 
classes,  conifer  and  its  three  seedling  classes,  annuals,  peren- 
nials, shrubs,  and  number  of  mature  conifers.  Those  computed 
include  total  herbs  and  other  combinations  of  compartments 
with  assigned  initial  values. 

Computation  of  flows  requires  that  ratios  of  present  value  to 
maximum  possible  value  for  numerous  compartments  be  com- 
puted at  each  time  step  (code  24700-27700).  The  FORTRAN 
name  for  such  a  ratio  is  the  FORTRAN  name  for  the  compart- 
ment with  a  prefix  of  P;  for  instance,  PASP  stands  for  the 
ratio  of  aspen  present  to  aspen  possible.  Should  such  a  ratio 
exceed  1,  the  value  of  the  corresponding  inhibitor  function 
becomes  negative,  which  causes  anomalous  behavior;  so  a 
check  is  made  to  reset  the  ratio  to  1  in  case  its  computed  value 
exceeds  1. 
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Regeneration  and  mortality  of  ASP1  (code  28000-29600)  are 
diagramed  in  figure  2.  ASP1  mortality  is  given  by  FLOW  15  = 
ASP1  *  P(41),  where  P(41)  is  the  site-specific  mortality  rate. 
ASP1  regeneration  is  given  by  FLOW  14  =  P(40)  *  Zl  *  Z2  * 
Z3  *  Z4  *  Z5  *  Z6  *  Z7  where  P(40)  is  the  site  specific  rate  and 
the  factors  Zl,  Z2,  Z3,  Z4,  Z5,  Z6,  and  Z7  are  the  restrictions 
contributed  by  ASP,  CON,  SHU,  ASP1,  ASP2,  ASP3,  and 
CONN  respectively  (fig.  3).  Note  that  as  the  ratio  of  present  to 
maximum  possible  biomass  goes  from  0  to  1 ,  the  corresponding 
restriction  increases  from  1  (no  restriction)  to  a  fraction 
nearer  0. 
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Figure  2.— Flow  diagram  depicting  regeneration  and  mortality  of  ASP1. 
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PASP  3 

FLOW  14  =  P(40)   *  Zl  *  Z2  *  Z3  *  Z4  *  Z5  *  Z6  *  Z7 

Figure  3.— Functions  that  inhibit  the  regeneration  of  ASP1. 
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Each  year  some  fraction  of  the  ASP1  is  allowed  to  graduate 
to  ASP2  and  a  fixed  proportion  of  the  ASP2  die  (code 
29700-31300).  The  flows  are  diagramed  (fig.  4)  and  the  nature 
of  the  restrictions  is  shown  in  figure  5.  The  ASP1  restriction  is 
considered  negligible  in  this  case,  but  the  herbs  are  not.  While 
the  herbs  may  not  restrict  suckering  of  aspen,  they  do  compete 
with  the  young  suckers  for  light. 
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Figure  4.—  Flow  diagram  depicting  graduation  from  ASP1  to  ASP2  and 
mortality  from  ASP2. 
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FLOW  16=  ASP1  *  P(51)  *  Zl  *  Z2  *  Z3  *  Z4  *     Z5  *  Z6  *  Z7 


Figure  5.— Functions  that  inhibit  the  graduation  of  ASP1. 


Graduation  from  ASP2  to  ASP3  and  mortality  of  ASP3 
(code  31400-32800)  are  shown  in  figure  6.  At  this  point  restric- 
tions due  to  herbs  and  ASP2  are  considered  negligible.  Inhibi- 
tions from  the  other  compartments  are  graphed  in  figure  7. 


Figure  6.—  Flow  diagram  depicting  graduation  from  ASP2  to  ASP3 
and  mortality  from  ASP3. 
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PASP  3 

FLOW  18  =  ASP2  *  P(59)   *  Zl  *  Z2  *  Z3  *  Z4  *  Z5 
Figure  7.— Functions  that  inhibit  the  graduation  of  ASP2. 
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To  ensure  a  conservative  system,  graduation  from  ASP3 
(code  32900-33600)  is  shown  in  figure  8  as  going  to  the  sink. 
The  sum  of  the  values  in  the  three  aspen  sucker  compartments 
and  their  source  and  sink  should  remain  constant.  Only  the 
conifer,  aspen,  and  shrub  compartments  have  identifiable 
restrictions  on  the  flow  (fig.  9). 
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Figure  8.— Flow  diagram  depicting  graduation 
of  ASP3. 
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FLOW  20  =  ASP 3  *  P(65)  *  Zl  *  Z2  *  Z3 

Figure  9.— Functions  that  inhibit  the  graduation  of  ASP3. 
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There  are  two  flows  from  the  source  to  ASP  in  figure  10. 
The  lower  flow  is  simply  the  conversion  of  ASP3  graduation 
from  numbers  to  biomass  (code  33700-34100).  The  upper  flow 
from  the  source  and  the  flow  to  the  sink  represent  aspen 
growth  and  mortality  (code  34200-35300).  The  restrictions  of 
the  aspen  and  conifer  on  aspen  growth  are  shown  in  figure  1 1 . 
The  restriction  due  to  aspen  itself  is  such  that  the  product  of 
the  restriction  with  the  growth  rate  would  equal  the  mortality 
rate  slightly  before  maximum  biomass  is  reached.  This  produces 
a  net  mortality  as  aspen  biomass  approaches  maximum  so  that 
aspen  biomass  will  not  continue  to  increase  past  maximum  due 
to  graduation  from  ASP3. 


CON 


Figure  10.— Flow  diagram  depicting  regeneration,  growth, 
and  mortality  of  aspen  (ASP). 
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Regeneration  and  mortality  of  CON1  (code  35400-37600)  are 
diagramed  in  figure  12.  As  in  the  case  of  ASP1 ,  the  relative 
mortality  and  maximum  rate  of  regeneration  are  site  specific. 
All  compartments  except  CON1  affect  the  regeneration  rate, 
though  annuals  and  perennials  are  pooled  as  a  single  effect  in 
HERB  and  the  aspen  suckers  are  converted  to  biomass  and 
pooled  with  shrubs  as  SHUH  (fig.  13).  The  restriction  of  aspen 
on  CON1  regeneration  is  the  only  inhibiting  function  which  is 
not  strictly  decreasing.  This  deviation  from  the  norm  was 
designed  to  account  for  a  "nurse"  effect  of  aspen  on  conifer 
regeneration.  Controversy  revolves  around  this  point,  however, 
with  some  arguing  that  there  is  no  positive  effect  from  aspen 
on  the  conifer  seedlings. 
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Figure  12.— Flow  diagram  depicting  regeneration  and  mortality 
ofCONL 
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Figure  13.— Functions  that  inhibit  the  regeneration  of  CON1. 
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Graduation  of  CONl  to  CON2  and  mortality  of  CON2 
(code  37500-39000)  are  diagramed  in  figure  14.  The  graduation 
is  affected  by  the  same  compartments  as  is  regeneration 
(fig.  15). 
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Figure  14.— Flow  diagram  depicting  graduation  of  CON1  to  CON2 
and  mortality  of  CON2. 
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FLOW  23  =  CON1  *  P(89)   *  Zl  *  Z2  *  Z3  *  Z4  *    Z5  *  Z6 
Figure  15.— Functions  that  inhibit  the  graduation  of  CON1. 
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The  diagram  of  CON2  graduation  and  CON3  mortality  (code 
3910O-4MO0)  is  in  figure  16.  In  the  graduation  of  CON2  the 
herbs,  shrubs,  and  aspen  suckers  are  no  longer  considered 
significant  factors.  The  remaining  inhibitions  are  shown  in 
figure  17. 


ASP 


Figure  16.— Flow  diagram  depicting  graduation  of 
CON2  to  CON3  and  mortality  of  CON3. 


PCON  2  PCON  3 

FLOW  25  =  C0N2  *  P(9  7)  *  Zl  *  Z2  *  Z3  *  Z4 
Figure  17.— Functions  that  inhibit  the  graduation  of  CON2. 


Like  ASP3,  the  C0N3  seedlings  graduate  to  a  numbers  sink 
(code  40500-41 100)  (fig.  18)  with  restrictions  (fig.  19).  The  sum 
of  the  values  in  the  source,  the  sink,  and  the  three  seedling 
compartments  should  remain  constant  over  the  years.  Because 
the  mortality  applied  to  the  conifer  numbers  (code 
42400-42800)  is  the  same  fraction  as  that  applied  to  the  conifer 
biomass,  all  the  loss  of  smaller  limbs  shows  up  as  loss  of  tree 
numbers,  so  that  the  numbers  estimate  grows  continually 
smaller  and  average  tree  size  estimates  are  too  large. 
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Figure  18.— Flow  diagram  depicting 
graduation  of  CON3. 
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Figure  19.— Functions  that  inhibit  the  graduation  of  CON3. 
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The  two  flows  into  the  conifer  compartment  from  the  source 
and  the  single  flow  out  to  the  sink  are  shown  in  figure  20.  The 
lower  flow  from  the  source  is  the  conversion  of  CON3  gradua- 
tion to  biomass  (code  41200-41600).  The  outflow  is  mortality 
(code  42900-43300),  and  the  upper  flow  is  in  conifer  growth 
(code  41700-42300),  which  is  affected  only  by  aspen  and  the 
conifer  itself  (fig.  21).  The  manner  in  which  aspen  affects  con- 
ifer growth  is  quite  critical  and  parameter  103,  which  gives  the 
maximum  restriction  of  aspen  on  conifer  growth,  is  probably 
the  one  to  which  the  model  is  most  sensitive.  Changing  this 
curve  can  drastically  alter  the  response  of  the  system.  The  in- 
hibition of  CON  on  CON  growth  is  designed,  as  in  the  case 
with  aspen,  to  prevent  an  increase  beyond  maximum  biomass  in 
the  CON  compartment  due  to  influx  from  CON3  graduation. 
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Figure  20.— Flow  diagram  depicting  regeneration,  growth,  and 
mortality  of  conifer  (CON). 
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FLOW  5  =  CON  *  P(106)   *  Zl  *  Z2 
Figure  21  .—Functions  that  inhibit  the  growth  of  conifer. 
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Production  and  mortality  of  annuals  (code  43400-44700)  are 
diagramed  in  figure  22.  All  other  compartments  affect  produc- 
tion, but  the  conifer  seedling  effects  are  pooled  into  CONN 
and  the  aspen  suckers  act  together  with  the  shrubs.  These 
restrictions  are  graphed  in  figure  23.  Because  the  time  step  is  1 
year,  the  relative  mortality  is  1 .  Sensitivity  of  the  system  to  the 
various  parameters  in  this  portion  of  the  program  is  very  slight, 
and  the  herbs  should  probably  not  be  subdivided  into  the  two 
compartments. 
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Figure  22.— Flow  diagram  depicting  production  and  mortality 
of  annuals. 
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Figure  23.— Functions  that  inhibit  the  production  of  annual  herbs. 
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Perennial  production  and  mortality  flows  (code  44800-46100) 
are  shown  in  figure  24.  Since  only  the  aboveground  biomass  is 
considered,  mortality  is  considered  to  be  100  percent  per  an- 
num. Restrictions  on  production  are  graphed  in  figure  25. 
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FLOW  9  =  P(122)   *  Zl  *  Z2  *  Z3  *  Z4  *  Z5 
Figure  25.— Functions  that  inhibit  the  production  of  perennial  herbs. 
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Regeneration,  growth,  and  mortality  (code  46200-48800)  of 
shrubs  are  diagramed  in  figure  26  with  restrictions  shown  in 
figures  27  and  28.  Regeneration  is  on  a  biomass  basis,  rather 
than  a  numbers  basis  as  it  was  for  aspen  and  conifer  reproduc- 
tion. All  compartments  including  the  shrub  compartment  affect 
the  regeneration  rate.  All  compartments  also  inhibit  the  growth 
of  shrubs,  though  effects  of  some  are  pooled.  The  restriction  of 
shrubs  on  shrub  growth  is  similar  to  that  of  aspen  and  conifer 
on  their  respective  growths — net  growth  is  blocked  just  prior  to 
maximum  possible  biomass  so  that  new  input  from  regenera- 
tion will  not  raise  the  biomass  beyond  the  maximum  stated,  as 
might  happen  if  aspen  and  conifer  were  absent  from  the  site. 

After  all  the  flows  for  one  time  step  have  been  computed, 
time  and  the  compartment  values  are  updated  (code 
48900-53200).  Thus,  all  values  in  one  step  are  completely  deter- 
mined by  the  values  in  the  preceding  step,  except  that  time  may 
influence  CON1  regeneration  because  a  few  years  are  required 
for  seedlings  to  become  established. 
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FLOW  11  =  P(130)   *  Zl  *  Z2  *  Z3  *  Z4  *  Z5  *  Z6 
Figure  27.— Functions  that  inhibit  the  regeneration  of  shrubs. 
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FLOW  12  =  SHU  *  P(135)   *  Zl  *  Z2  *  Z3  *  Z4  *  Z5  *  Z6 
Figure  28.— Functions  that  inhibit  the  growth  of  shrubs. 
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MODEL  RESULTS 
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Figure  29a.— Biomass  of  four  major 
system  components  graphed  against 
time  for  "standard  run"  in  which  all 
parameters  are  at  nominal  levels. 
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Figure  29b.— Density  of  three  classes 
of  aspen  suckers  graphed  against 
time  for  "standard  run"  in  which  all 
parameters  are  at  nominal  levels. 


The  model  results  are  quite  realistic.  Although  no  one  has 
directly  observed  400  years  of  succession  in  the  Intermountain 
West,  experts  agree  that  the  model  output  for  a  typical  site 
fits  their  expectations  developed  from  observations  of  a  series  of 
existing  sites. 

Results  of  the  model  for  400  years  following  a  hot  fire  are 
shown  in  figure  29  (a,  b,  c).  In  this  "standard  run  "  it  is 
assumed  that  maximum  biomass  possible  for  the  major 
compartments  is  200  000  kg/ha  for  aspen,  250  000  kg/ha  for 
conifer,  10  000  kg/ha  for  shrubs,  and  4  000  kg/ha  for  herbs. 
While  peak  aspen  biomass  (191  000  kg/ha)  occurs  in  the  132d 
year,  the  aspen  biomass  exceeds  190  000  kg/ha  for  approx- 
imately 40  years.  The  conifers  come  in  much  more  slowly  and 
conifer  biomass  does  not  exceed  aspen  biomass  until  nearly  300 
years  after  the  disturbance.  After  400  years  the  conifer  biomass 
is  202  000  kg/ha  and  still  increasing  slightly.  The  shrubs  peak 
in  the  80th  year  with  7  050  kg/ha  and  herbs  peak  the  eighth 
year  with  a  biomass  of  2  700  kg/ha.  Graphs  of  biomass  versus 
time  for  these  major  compartments  are  shown  in  figure  29a. 
Aspen  suckers  come  in  quite  rapidly  (fig.  29b),  peaking  within 
the  first  few  years,  but  then  they  decline  rapidly.  The  smallest 
size  class  (ASP1)  is  the  most  numerous.  The  number  of  conifer 
seedlings  gradually  increases  to  a  peak  value  about  300  years 
following  the  disturbance  and  then  decreases  again,  the  small 
seedlings  being  more  numerous  (fig.  29c). 
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Figure  29c— Density  of  three  classes 
of  conifer  seedlings  graphed  against 
time  for  "standard  run"  in  which  all 
parameters  are  at  nominal  levels. 
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By  changing  the  maximum  rate  at  which  conifer  seedlings 
become  established  (parameter  79),  from  1,000  to  10,000,  an 
earlier  dominance  of  conifers  is  obtained  (fig.  30).  This  results 
in  an  earlier  suppression  of  the  other  compartments  and  a 
shorter  time  period  for  conifers  to  approach  maximum. 
Although  such  behavior  may  occur  in  the  Intermountain  area, 
it  is  probably  more  typical  of  aspen  woodlands  of  the  Upper 
Midwest. 
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Figure  30.— Biomass  of  four  major 
system  components  graphed  against 
time  for  case  simulating  rapid 
regeneration  of  conifers. 


If  other  disturbances  are  not  allowed  and  if  conifer  seedlings 
do  not  become  established  (no  seed  source),  the  values  in  the 
other  compartments  remain  much  higher  (fig.  31)  throughout 
the  simulation  run. 

In  figure  32  an  aspen  harvest  is  simulated  at  90  years  in  the 
successional  process  to  favor  conifer  production  on  the  site. 
Conifer  growth  is  shown  to  be  enhanced  by  that  harvest;  in 
just  270  years  it  attains  the  value  achieved  in  400  years  in  the 
"standard  run,  "  and  after  400  years  it  is  close  to  its  maximum 
biomass.  Immediately  after  the  aspen  harvest,  the  shrubs  and 
herbs  increase  but  are  quickly  suppressed  again  by  the  rapidly 
growing  conifers.  Aspen  comes  in  again  after  the  harvest  but 
peaks  at  only  85  000  kg/ha  and  is  almost  eliminated  by  the 
400th  year.  Aspen  regeneration  is  stimulated  by  the  harvest  of 
aspen  overstory,  but  establishment  of  conifer  seedlings  is  sup- 
pressed (fig.  33).  This  situation  brings  up  a  controversial  point 
in  the  model:  the  effect  of  aspen  on  conifer  seedling  regenera- 
tion. Foresters  differ  as  to  the  influence  of  aspen  on  conifer 
regeneration,  and  it  was  decided  to  use  a  positive  influence  in 
this  model.  This  positive  effect  on  conifer  regeneration  is  par- 
tially due  to  the  ameliorating  effects  aspen  has  on  soil  surface 
temperatures  (Sperger  1980).  It  is  accomplished  by  having  the 
aspen  restriction  on  conifer  regeneration  decrease  as  aspen 
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Figure  31.— Biomass  of  four  major 
system  components  graphed  against 
time  for  case  in  which  there  is  no 
conifer  regeneration. 
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Figure  32.— Biomass  of  four  major 
system  components  graphed  against 
time  for  case  in  which  aspen  is 
harvested  at  90  years. 
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Figure  33.— Density  of  conifer  seedlings 
in  case  where  aspen  overstory  is 
harvested  at  90  years. 

biomass  increases.  Thus,  harvest  of  the  aspen  overstory  imposes 
a  restriction  on  the  regeneration  of  conifer.  This  influence 
merits  further  investigation.  However,  even  if  no  influence  of 
aspen  on  conifer  regeneration  is  allowed,  the  conifer  biomass  in 
the  300th  year  is  only  1.9  percent  higher  for  the  "standard  run," 
while  peak  biomass  for  shrubs,  aspen,  and  herbs  changes  very 
little. 

A  more  typical  clearcutting  situation  is  shown  in  figure  34. 
Here  all  the  overstory  conifers  and  half  the  conifer  seedlings 
are  removed  at  300  years.  Also,  half  the  mature  aspen  biomass 
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Figure  34.— Biomass  of  four  major 
system  components  graphed  against 
time  for  case  in  which  conifer  and 
aspen  are  harvested  at  300  years. 


is  removed  by  harvesting.  The  next  300  years  are  very  similar  to 
the  first  300  except  that  aspen  peaks  a  little  more  quickly  and 
at  a  somewhat  higher  value.  One  should  bear  in  mind  that  this 
second  peak  does  not  represent  the  even-aged  stand  represented 
by  the  first  peak. 

The  major  shortcoming  of  the  model  is  lack  of  confidence  in 
the  parameter  values  used.  It  is  quite  possible  that  widely  dif- 
fering sets  of  parameters  can  give  very  similar  model  results. 
On  the  other  hand,  slightly  differing  parameters  can  yield  wide- 
ly differing  results.  To  demonstrate  the  degree  of  uncertainty  in 
the  model,  32  simulation  runs  (variations  on  the  "  standard  run" ) 
were  made.  Parameters  for  the  maximum  production  rates  of 
annuals  and  perennials  and  for  the  maximum  values  of  the 
aspen,  conifer,  shrub,  annual,  and  perennial  compartments 
were  fixed  for  the  "typical"  site  investigated  and  the  initial 
value  for  number  of  conifers  was  held  at  zero.  At  the  start  of 
each  of  these  32  runs,  each  of  the  remaining  130  parameter 
values  was  picked  at  random  from  a  uniform  distribution  over 
the  range  of  values  considered  possible  for  that  parameter  (ap- 
pendix C)  and  held  for  the  duration  of  the  run.  Choice  of 
parameters  was  not  always  completely  random,  however; 
checks  were  made  to  ascertain  that  the  sum  of  relative  flow 
rates  from  a  compartment  did  not  exceed  1  and  that  aspen 
growth  rate  never  exceeded  the  conifer  growth  rate.  The  basic 
shapes  of  the  resulting  graphs  were  much  the  same  in  all  cases. 
That  is,  if  appropriately  different  scales  were  used  for  the  32 
graphs  of  any  major  compartment,  they  would  be  quite  similar 
in  appearance.  Some  idea  of  the  distribution  of  these  32  runs 
can  be  gained  from  figure  35,  which  shows  for  each  point  in 
time  the  minimum  and  maximum  values  of  each  of  the  four 
major  compartments  as  well  as  the  25th,  50th,  and  75th  percen- 
tile values. 
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Figure  35.— Results  from  32  com- 
puter runs  with  randomly  varying 
parameters  showing  maximum,  75th 
percentile,  50th  percentile,  25th 
percentile,  and  minimum  values  at 
each  point  in  time. 
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Sensitivity  Analysis 

To  test  parameter  sensitivity,  the  same  130  parameters  men- 
tioned in  the  preceding  section  were  altered  one  at  a  time  from 
the  nominal  value  to  the  more  remote  end  point  in  the  range  of 
possible  values.  Change  in  the  peak  value  for  each  of  the  four 
major  compartments  was  recorded,  as  was  the  change  in  time 
at  which  the  peak  occurred.  This  exercise  showed  that  the 
model  is  fairly  insensitive  to  changes  in  many  of  the  para- 
meters, indicating  that  the  model  may  be  overdetermined;  how- 
ever, it  is  still  quite  sensitive  to  changes  in  others.  The  change 
in  time  at  which  peaks  occurred  varied  little,  and  the  peak 
biomass  of  herbs  also  did  not  indicate  much  sensitivity  except 
to  the  initial  value  of  annuals.  Parameters  to  which  the  peak 
values  of  aspen,  conifers,  and  shrubs  are  most  sensitive  are 
those  involving  the  growth  and  mortality  rates  of  those  com- 
partments, the  restriction  of  aspen  on  conifer  growth,  the 
restriction  of  conifers  on  aspen  growth,  and  regeneration  into 
or  mortality  from  CON1 .  For  instance,  an  increase  in  the  an- 
nual growth  rate  of  conifers  from  7  percent  to  7.5  percent 
reduces  the  peak  value  of  aspen  biomass  by  0.7  percent  and  in- 
creases the  peak  conifer  biomass  by  19.7  percent.  If  the  max- 
imum restriction  of  aspen  on  conifer  growth  is  lessened  from 
0.5  to  0.7,  the  peak  aspen  biomass  is  decreased  3.1  percent 
while  the  peak  conifer  biomass  is  increased  70.9  percent. 


Further  investigations  of  these  relationships  indicated  that 
one  could  raise  conifer  growth  rates  from  7  percent  to  8  per- 
cent, raise  conifer  mortality  from  2  to  3  percent,  lighten  the 
restriction  of  aspen  on  conifer  growth  from  0.5  to  0.5625,  and 
get  almost  no  change  at  all  in  the  model  results!  Therefore,  the 
parameter  values  are  not  unique,  but  it  appears  that,  so  long  as 
there  is  some  method  of  bringing  growth  and  mortality  rates 
together,  the  conifer  biomass  will  show  appropriate  sigmoidal 
growth. 

Some  two-way  interactions  were  also  tested.  The  three  in- 
dicators chosen  were  peak  aspen  and  peak  shrub  biomass  in 
400  years  and  conifer  biomass  at  300  years.  Because  graphs  of 
a  compartment's  size  are  very  similar  if  scales  are  adjusted  ap- 
propriately, it  was  felt  that  maximum  values  would  be  good  in- 
dicators. Because  the  peak  of  herbs  is  almost  never  affected,  it 
was  dropped  as  an  indicator.  Peak  conifer  biomass  always  oc- 
curred at  400  years  and  differences  at  that  time  could  be  slight 
while  differences  earlier  were  much  greater;  so  biomass  at  300 
years  was  a  more  appropriate  indicator  for  conifer  biomass. 

Of  the  130  parameters  allowed  to  vary  it  was  found  that  32 
caused  at  least  1  of  those  3  indicators  to  change  from  its 
nominal  value  by  at  least  2  percent  (table  1).  For  these  32 
parameters,  all  496  paired  interactions  were  evaluated.  For  I, 
1(A),  1(B),  and  I(AB)  (defined,  respectively,  as  the  values  of  the 


Table  1.— Parameters  (identified  in  appendix  C)  to  which  the  ASPEN  model  is  most  sensitive,  their  nominal  and 
perturbed  values,  and  the  effect  of  those  perturbations  on  the  three  indicators 


Parameter 

Nominal 

Perturbed 

Peak  aspen 

Conifer  at  300 

Peak  shrub 

number 

value 

value 

%  change 

%  change 

%  change 

13 

15000. 

20000. 

-0.3 

2.8 

-2.0 

67 

.03 

.035 

-1.9 

4.7 

2.0 

70 

1.2 

1.4 

1.2 

-2.6 

-.1 

77 

.2 

.4 

-.4 

3.6 

-.4 

78 

.001 

.2 

-.7 

3.8 

-.8 

79 

1000. 

1250. 

-.9 

7.6 

-.8 

80 

5. 

10. 

.4 

-3.4 

.7 

81 

.7 

.9 

.6 

-7.0 

.6 

82 

.5 

.7 

-.3 

4.1 

-.2 

87 

.7 

.9 

-.2 

2.3 

-.2 

88 

.7 

.9 

-.3 

2.5 

-.3 

89 

.3 

.4 

-.9 

7.8 

-.8 

90 

.45 

.6 

.6 

-7.1 

.6 

91 

.5 

.7 

-.3 

3.8 

-.2 

97 

.25 

.4 

-1.3 

10.8 

-1.2 

98 

.1 

.15 

.5 

-6.0 

.5 

99 

.5 

.7 

-.2 

2.6 

-.2 

101 

.2 

.3 

-.6 

5.8 

-.9 

102 

6. 

7. 

-.6 

3.7 

-.7 

103 

.5 

.7 

-3.1 

70.9 

-1.5 

104 

3. 

4. 

.4 

-13.6 

.3 

105 

.02 

.025 

1.0 

-42.5 

.6 

106 

.07 

.075 

.7 

19.7 

.4 

115 

.5 

.7 

.3 

-2.9 

-4.8 

119 

.5 

.4 

3.1 

-3.1 

-2.6 

123 

4. 

6. 

.0 

.4 

-2.3 

130 

100. 

125. 

.1 

-.5 

2.6 

131 

.5 

.6 

.0 

-.8 

8.3 

133 

.5 

1. 

.1 

-1.5 

10.6 

134 

.03 

.04 

-.1 

1.9 

-17.9 

135 

.2 

.25 

.2 

-2.0 

14.2 

138 

.4 

.6 

.1 

-1.2 

9.3 

30 


indicator  with  no  parameters  perturbed,  only  parameter  A  per- 
turbed, only  parameter  B  perturbed,  and  both  A  and  B  per- 
turbed), the  effect  of  the  interaction  of  A  with  B  is  defined  to 
be  [I(AB)  - 1]  -  [1(A)  -  I]  -  [1(B)  -  I].  In  only  24  of  these  496 
cases  did  effect  of  the  interaction  on  an  indicator  exceed  2  per- 
cent of  the  nominal  value  of  the  indicator.  In  every  case  it  was 
the  300-year  biomass  of  the  conifer  that  was  changed  by  at 
least  2  percent.  In  20  of  the  24  cases  the  parameter  giving  the 
maximum  restriction  of  aspen  on  conifer  growth  was  involved. 
The  other  four  cases  involved  growth  or  mortality  of  aspen  or 
conifer.  In  none  of  the  cases  investigated  did  the  effect  of  the 
interaction  exceed  the  larger  of  the  two  main  effects.  Because 
parameters  having  the  greatest  effect  on  the  indicators  were  in- 
vestigated, it  is  likely  that  a  similar  result  would  have  held  if 
the  other  5,000+  pairs  had  been  checked.  Thus,  the  strongest 
interactions  affecting  the  selected  indicators  probably  are  the  24 
listed  in  table  2. 

The  main  effects  and  interactions  previously  mentioned  in- 
volved altering  just  one  or  two  parameters  at  a  time.  Another 
concept  of  effect  involves  a  factorial  design  in  which  all  com- 
binations of  parameters  are  perturbed  and  the  average  effect 
over  all  combinations  of  other  factors  is  considered.  (A  discus- 
sion of  factorial  designs  and  fractional  factorial  designs  can  be 
found  in  Cochran  and  Cox,  1975.)  A  complete  factorial  design 
for  130  factors,  or  even  30  factors,  would  have  been  too  expen- 
sive; however,  a  fractional  factorial  design  for  a  set  of  36 
parameters,  which  included  the  32  previously  found  to  have  the 
largest  main  effects,  was  possible.  These  were  divided  into  six 


groups  called  macroparameters  because  each  is  treated  as  a  unit 
in  the  factorial  analysis.  These  macroparameters,  named  A,  B, 
C,  D,  E,  and  F,  contain  six  parameters  each.  Care  was  taken 
to  separate  closely  related  parameters  into  different  macro- 
parameters.  Then  a  0.5  replication  of  a  2**6  factorial  design 
was  done  using  ABCDEF  as  the  defining  contrast.  The  three 
indicators  of  interest  were  again  the  peak  values  of  aspen  and 
shrubs  over  400  years  and  the  value  of  conifer  at  300  years. 
Contents  of  the  macroparameters  are  shown  in  table  3.  The 
ones  perturbed  in  each  repetition  of  the  experiment  and  cor- 
responding values  of  the  three  indicators  are  shown  in  table  4. 

The  effect  of  a  factor  on  an  indicator  is  defined  to  be  1/32 
of  the  sum  of  all  values  of  that  indicator  obtained  from  experi- 
ments in  which  an  odd  number  of  the  macroparameters  within 
that  factor  were  perturbed,  minus  1/32  of  the  sum  of  all  values 
of  the  indicator  obtained  from  experiments  in  which  an  even 
number  of  the  macroparameters  within  that  factor  were  per- 
turbed. Each  main  effect  has  an  alias  with  a  5-way  interaction; 
each  2-way  interaction  has  an  alias  with  a  4-way  interaction; 
and  each  3-way  interaction  has  an  alias  with  another  3-way  in- 
teraction. The  various  factors,  their  aliases,  corresponding 
change,  and  percentage  change  in  values  of  the  three  indicators 
are  given  in  table  5.  From  this  table  one  can  observe  that  for 
all  three  indicators  the  effects  of  the  1-way  (5-way)  interactions 
are  in  general  greater  than  the  effects  of  the  2-way  (4-way)  in- 
teractions, which  are  in  general  larger  than  the  effects  of  the 
3-way  interactions.  This  observation  and  intuition  indicate  that 
the  larger  effects  are  associated  with  the  lower  rather  than  the 


Table 2.— Parameter  pairs  to  which  the  system  is  most  sensitive,  their  nominal  and  perturbed  values,  and  the 
effect  of  each  perturbation,  due  to  the  interaction,  on  the  conifer  biomass  in  the  300th  year 


Percent  change 
in 

Parameter         Nominal         Perturbed        Parameter        Nominal        Perturbed  CON  (300) 

number  value  value  number  value  value  due 

to  interaction 


13 

15000. 

20000. 

103 

0.5 

0.7 

-2.2 

67 

.03 

.035 

104 

3. 

4. 

-2.2 

67 

.03 

.035 

105 

.02 

.025 

-2.4 

67 

.03 

.035 

106 

.07 

.075 

2.1 

77 

.2 

.4 

103 

.5 

.7 

-2.7 

78 

.001 

.2 

103 

.5 

.7 

-3.7 

79 

1000. 

1250. 

103 

.5 

.7 

-5.6 

80 

5. 

10. 

103 

.5 

.7 

2.1 

81 

.7 

.9 

103 

.5 

.7 

5.1 

82 

.5 

.7 

103 

.5 

.7 

-3.3 

89 

.3 

.4 

103 

.5 

.7 

-5.8 

90 

.45 

.6 

103 

.5 

.7 

5.2 

91 

.5 

.7 

103 

.5 

.7 

-3.1 

97 

.25 

.4 

103 

.5 

.7 

-8.1 

98 

.1 

.15 

103 

.5 

.7 

4.5 

99 

.5 

.7 

103 

.5 

.7 

-2.1 

101 

.2 

.3 

103 

.5 

.7 

-4.3 

102 

6. 

7. 

103 

.5 

.7 

-3.8 

103 

.5 

.7 

104 

3. 

4. 

8.1 

103 

.5 

.7 

105 

.02 

.025 

23.3 

103 

.5 

.7 

106 

.07 

.075 

-14.5 

103 

.5 

.7 

119 

.5 

.7 

2.1 

103 

.5 

.7 

115 

.5 

.7 

2.1 

104 

3. 

4. 

105 

.02 

.025 

3.8 

31 


higher  order  interactions.  In  other  words  this  fractional  fac- 
torial design  serves  to  reinforce  the  previous  results  determining 
which  parameters  the  system  is  most  sensitive  to  without  giving 
as  much  detail.  The  effect  of  a  macroparameter  corresponds 
roughly  to  the  effects  of  its  constituent  parts  but  is  not  the 
mean  of  those  separate  effects. 


Table  3.— Elements  of  a  macroparameter  used  in  the  fractional 
factorial  design 


Table  4.— Treatment  combinations  in  the  fractional  factorial 

design  and  corresponding  values  of  the  three  indicators 


Macroparameter 


Parameters 


A 

36 

68 

97 

104 

115 

135 

B 

70 

80 

82 

98 

130 

134 

C 

13 

67 

77 

88 

99 

123 

D 

58 

78 

89 

101 

106 

133 

E 

32 

79 

90 

102 

103 

131 

F 

71 

81 

91 

105 

119 

138 

Conifer 

Peak 

biomass 

Peak 

Treatment 

aspen 

at  300 

shrub 

combinations 

biomass 

years 

biomass 

none 

191431. 

136588. 

7050. 

AB 

195503. 

118820. 

6928. 

AC 

^  O 7AOO 

187023. 

143652. 

7497. 

AD 

187145. 

165280. 

8307. 

AE 

185432. 

230285. 

8005. 

AF 

195483. 

62685. 

8181. 

BC 

189827. 

147357. 

5792. 

BD 

190058. 

172190. 

6975. 

DC 

188938. 

230019. 

6646. 

br 

19b33o. 

67952. 

6944. 

UD 

172764. 

21 1782. 

7859. 

CE 

171980. 

243178. 

7198. 

Or 

18922b. 

92279. 

7522. 

Ub 

^  "7-1  noc 

1  f  iy2b. 

244/91. 

8316. 

L)r 

-i  on  a  on 

121  lid. 

8175. 

br 

188281. 

206687. 

7974. 

ABCD 

184586. 

175282. 

7593. 

AbUb 

183182. 

235315. 

7121. 

AbOr 

iy4Jb9. 

80309. 

7514. 

AbUb 

18o  loo. 

2>j8d80. 

8150. 

ABDF 

194563. 

105379. 

8065. 

ABEF 

193766. 

190216. 

7968. 

AOUb 

1  CO  1 oc 

loo  180. 

247128. 

8679. 

AGUr 

182yoo. 

12y3/3. 

8621. 

AObr 

181381. 

210800. 

8325. 

ADEF 

181111. 

225440. 

8884. 

BCDE 

167716. 

246617. 

7526. 

BCDF 

186895. 

134833. 

7472. 

BCEF 

185712. 

214236. 

7136. 

BDEF 

185898. 

225754. 

7990. 

CDEF 

162565. 

243121. 

8545. 

ABCDEF 

177484. 

232844. 

8444. 

Table  5.— Factors,  their  aliases,  and  interactive  effects  on  the  three  indicators  as  determined  by  the  fractional 
factorial  design 

Defining 

contrast  Peak  aspen  Conifer  biomass  Peak  shrub 

ABCDEF  biomass  at  300  years  biomass 


%  from  %  from  %  from 

Factor  Alias  Effect         nominal  Effect         nominal  Effect  nominal 


A 

BCDEF 

1135. 

0.6 

-4458. 

-3.3 

286. 

4.1 

B 

ACDEF 

3177. 

1.7 

-3251. 

-2.4 

-340. 

-4.8 

C 

ABDEF 

-4457. 

-2.3 

7830. 

5.7 

-54. 

-.8 

D 

ABCEF 

-4419. 

-2.3 

15778. 

11.6 

369. 

5.2 

E 

ABCDF 

-5028. 

-2.6 

50143. 

36.7 

200. 

2.8 

F 

ABCDE 

2394. 

1.3 

-19946. 

-14.6 

254. 

3.6 

AB 

CDEF 

-183. 

-.1 

576. 

.4 

45. 

.6 

AC 

BDEF 

270. 

.1 

-148. 

-.1 

10. 

.1 

AD 

BCEF 

235. 

.1 

-633. 

-.5 

-44. 

-.6 

AE 

BCDF 

157. 

.1 

2040. 

1.5 

-21. 

-.3 

AF 

BCDE 

-88. 

0. 

422. 

.3 

-21. 

-.3 

BC 

ADEF 

804. 

.4 

-469. 

-.3 

-13. 

-.2 

BD 

ACEF 

840. 

.4 

-318. 

-.2 

17. 

.2 

BE 

ACDF 

885. 

.5 

579. 

.4 

31. 

.4 

BF 

ACDE 

-391. 

-.2 

399. 

.3 

46. 

.7 

CD 

ABEF 

-1178. 

-.6 

-224. 

-.2 

46. 

.7 

CE 

ABDF 

-1187. 

-.6 

-2432. 

-1.8 

-6. 

-.1 

CF 

ABDE 

439. 

.2 

727. 

.5 

16. 

.2 

DE 

ABCF 

-1245. 

-.7 

-7113. 

-5.2 

16. 

.2 

DF 

ABCE 

442. 

.2 

2244. 

1.6 

-79. 

-1.1 

EF 

ABCD 

461. 

.2 

9827. 

7.2 

-27. 

-.4 

ABC 

DEF 

-38. 

0. 

1618. 

1.2 

2. 

0. 

ABD 

CEF 

-34. 

0. 

1113. 

.8 

-2. 

0. 

ABE 

CDF 

56. 

0. 

-604. 

-.4  . 

-12. 

-.2 

ABF 

CDE 

-196. 

-.1 

-793. 

-.6 

-4. 

-.1 

ACD 

BEF 

19. 

0. 

-1227. 

-.9 

-11. 

-.2 

ACE 

BDF 

-30. 

0. 

558. 

.4 

-5. 

-.1 

ACF 

BDE 

157. 

.1 

917. 

.7 

3. 

0. 

ADE 

BCF 

-48. 

0. 

1027. 

.8 

1. 

0. 

ADF 

BCE 

132. 

.1 

614. 

.4 

7. 

.1 

AEF 

BCD 

206. 

.1 

-1190. 

-.9 

3. 

0. 
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Use  of  the  Program 

To  run  the  program  ASPEN  one  must  specify  from  which 
device  to  read  the  data  and  on  which  device  to  print  the  results. 
Appendix  A  contains  a  listing  of  the  program  as  it  is  run  with  a 
card  deck  at  Utah  State  University  (USU)  and  the  input  device, 
file  5,  is  specified  as  the  card  reader  while  the  output  device  on 
file  6  is  specified  as  the  printer  (code  200-300). 

A  brief  flow  chart  of  the  program  is  given  in  figure  36.  It 
does  not  describe  the  flow  equations  in  detail  because  they  are 
handled  in  the  section  on  model  structure. 

The  first  information  required  from  the  data  file  concerns 
time.  The  values  of  DT,  TEND,  STEP,  and  PSTSZ  are  read 
and  then  printed  as  part  of  the  program  output  (code 
9300-9900).  DT  is  the  step  size  of  the  difference  equations. 
Although  it  is  a  parameter,  it  was  always  used  as  1  year  in  our 
work  at  USU.  While  it  is  conceivable  that  one  might  want  to 
increase  step  size  to  save  computation,  some  of  the  parameters 
would  need  to  be  altered  to  compensate.  For  instance,  a  growth 
rate  of  8  percent  per  year  is  not  quite  the  same  as  a  growth  rate 
of  16  percent  over  2  years,  and  the  difference  can  be  pro- 
nounced after  400  years.  The  length  of  the  run  in  years  is 
TEND;  the  number  of  years  between  printout  of  state  variables 
is  STEP;  and  PSTSZ  specifies  the  number  of  years  between 
points  on  the  graphical  output.  The  graphing  subroutines  (code 


56000-65200)  allow  for  only  101  points  (including  the  point  at 
time  0)  plotted  across  the  page;  so  PSTSZ  must  be  an  integer 
no  smaller  than  the  quotient  TEND  divided  by  100.  In  a  run  of 
100  years  or  fewer,  each  year's  values  can  be  graphed,  but  for 
a  run  of  400  years,  no  more  than  every  fourth  year  can  be 
graphed  (PSTSZ  =  400/100  =  4  or  larger).  The  printed  output 
is  fairly  voluminous,  and  printing  every  year's  values  can  be  ex- 
pensive. A  sample  output  is  shown  in  appendix  B. 

Immediately  following  the  time  information  are  the  138 
parameter  values.  The  program  reads  these  values  and  prints 
them  with  the  other  output  (code  10000-1 1000).  Parameter 
values  are  read  in  free  format,  five  to  the  line.  The  last  two 
values,  P(139)  =  P(140)  =  0,  are  simply  space  fillers.  It  is  not 
easy  to  find  a  particular  value  from  the  listing,  but  it  is  possible 
if  it  is  remembered  that  the  four  values  on  the  first  line  are  for 
time  and  that  the  parameter  values  follow  in  order,  five  to  the 
line. 

Time  is  always  initiated  at  zero  by  the  program.  If  the  pro- 
gram is  stopped  and  restarted  again  at  a  later  year,  the  later 
year  does  not  show  up  in  the  output;  time  is  reinitiated  to  zero. 
The  time  delay  for  conifer  seedling  establishment,  parameter 
80,  tells  how  many  years  into  the  run  to  delay.  For  example  if 
the  run  is  stopped  for  some  reason  and  restarted  at  100  years 
and  if  a  5-year  time  delay  on  conifer  seedling  establishment  is 
desired,  parameter  80  should  be  set  at  5,  not  105. 
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Figure  36.— A  brief  flow  chart  of  the  program  ASPEN. 
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CONCLUSIONS 

For  trajectory  analysis  of  orbiting  satellites,  a  flat-earth 
theory  is  woefully  inadequate,  but  for  consideration  of  baseball 
trajectories  at  a  ball  park,  that  theory  is  quite  adequate. 
Similarly,  the  succession  model  presented  in  this  paper,  while 
certainly  incomplete,  might  help  a  land  manager  identify  the 
successional  stage  of  a  forest  and  determine  how  its  composi- 
tion might  change  in  the  future. 

The  model  yields  believable  results  for  the  forest  system 
simulated,  whether  that  forest  was  undisturbed  or  subjected  to 
various  manipulations.  It  could  be  of  some  aid  in  making 
forest  management  decisions  if  managers  bear  in  mind  that  the 
model  is  a  partial  truth  and  not  the  whole  truth. 

The  distressing  shortcoming  of  the  model  is  that  good 
estimates  do  not  exist  for  most  of  the  parameters  used.  Many 
parameter  values  were  chosen  from  a  range  which  scientists  in 
our  project  estimated  to  include  the  actual  value.  To  some  of 
these  parameters  the  model  is  quite  sensitive,  giving  vastly  dif- 
ferent simulation  results  depending  on  which  value  in  that 
range  is  chosen.  One  should  not  rely  on  the  model  as  a  predic- 
tor until  closer  estimates  of  these  parameters  are  obtained. 

Sensitivity  analysis  has  also  shown  that  the  model  is  quite  in- 
sensitive to  many  parameter  estimates,  indicating  that  a 
simplification  effort  might  substantially  reduce  the  model  com- 
plexity while  keeping  the  essential  features.  We  have  also 
isolated  those  parameters  (growth  and  mortality  rates  of  aspen 
and  conifer  and  the  restrictive  effect  of  aspen  on  conifer 
growth)  that  have  the  most  drastic  effects  on  the  model  output. 
Further  research  in  these  areas  would  be  most  valuable  to  our 
understanding  of  the  system  and  would  lend  confidence  to  the 
model  as  a  predictive  tool. 
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C«,9ff«                          SUC  *  TOTAL  ASPEN  SUCKERS             #/HEC  OOOOSBon 

Cffnt                        SEED  a  TOTAL  CONIFER  SEEDLINGS     #/HEC  00005900 

CtftM                        CONN  a  MASS  OF  CONIFER  SEEDLINGS       KG /HE C  00006000 

C9>,99                    6IOSUH  a  SUH  OF  BIOMA5S  STATE  VARIABLES  00006100 

C,,...                                               CSHOULD  RE   ZERO )           KG/HEC  00006200 

Citttf                    SUCSUM  a  SUM  OF  ASPEN  SUCKER  STATE  VARIABLES  00006300 

C.att.                                               (SHOULD  P'E   ZERO)           #/HEC  00006400 

C...M                    SEESUM  a  SUM  OF  CONIFER  SEFDLING  STATE  VARIABLES  0000^500 

Ctctf                                               (SHOULD  BE   ZERO )           #VHEC  00006600 

Ctfiit                 AN  M  PRECEDING  ThE  VARIABLE  NAME  REPRESENTS  00006700 
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C««»tf  00006900 
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WRITE (6, 33) M ASP  1 #MASP2.MASP 3, MC ONi ,MCON2, ^CON3 
FORMAT ( 1 H   ,»MA3P1    a»fF10.2#2X#  'WASP2  a'EiO^^X*1 
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WRITE  (6.341MASPS.MPFR.  MANN 
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Ct... 
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C»  i  M 

c  •  •  • » 

c... 
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SET   TIME   AND  COUNT   TO  0.0 

TIME  a  o.O 
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ASP2  a  P(21) 
ASP3  s  P(22) 
CONi   a  P(23) 
C0N2  a  p(2a) 
CQN3  =  P(25) 
NOCQN  8  P(26) 


2) 


C  • . .  •  t 

c  

c  


UPDATE  CORRESPONDING  VALUES  OF  OTHER  VARIABLES 


HE«8 

sue 

SEED 
ASPS  = 

5HUH  a 
CONN  a 
BIOSOR 
8I0SIN 
SUCSOR 
SUCSIN 
SEESOR 
SEESIN 
BIOSUM 
SUCSUM 
SEE  SUM 


♦  P(29)*ASP3 


77 


ANN  ♦  PER 
ASpt*AsP2*AsP3 
C0NUC0N2^C0N3 
P(27)*ASP1   ♦  P(28)*ASP2 
SHU   ♦  ASPS 

PC30)*CON1    ♦  P(31)*C0N2  ♦  P(32)*C0N3 
a  «ASP  ■  CON  -   ANN   •  PER  ■  SHU 
a  0,0 

a  *ASPt   •  ASP2  -  ASP3 
a  0.0 

a  -CONI   m   C0N2  •  C0N3 
a  0.0 

a  BIOSOR  t  ASP  ♦  CON  t 
a  SUCSOR  ♦  ASPl   ♦  ASP2 
a  SEESOR  t  CONI   ♦  C0N2 
IF    (NOCqN.GT.0.0)   GO  TO  77 
AWPC   a  m\ 
GO  TO  78 

A*PC   a  CON/NOCOM 


ANN  ♦  PER  ♦  SHU 

♦  ASP3  +  SUCSIN 

♦  C0N3  ♦  SEESIN 


♦  8I0SIN 


0  00 

11900 

000 

L2000 

000 

12100 

0  00 

12200 

000  1 

12300 

000 

12400 

000" 

12500 

000 

12600 

000 

12700 

000' 

12800 

000 

12900 

000 

13000 

000 

13100 

000 

13200 

ooo 

13300 

000 

I  3400 

000 

13500 

ooo 

13600 

ooo 

13700 

ooo 

I  3800 

ooo 

1  3900 

ooo 

14000 

ooo 

14100 

000 

14200 

ooo 

L4300 

000 

14000 

000 

I  4500 

000 

14600 

ooo 

14700 

000 

14800 

ooo 

14900 

000 

15000 

000 

15100 

000 

15200 

000 

15300 

ooo 

15400 

000 

15500 

000 

15600 

ooo 

15700 

0001 

15800 

0001 

15900 

oool 

16000 

0001 

16100 

0001 

16200 

0001 

16300 

000  1 

16400 

0001 

16500 

oool 

16600 

oool 

16700 

0001 

6800 

oool 

-6900 

oool 

7000 

oool 

7100 

oool 

7200 

oool 

7300 

0001 

7400 

0001 

7500 

oool 

7600 

0O01 

7700 

0001 

7800 

0001 

7900 
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78  CONTINUE  00018000 

Cmmi  00018100 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCrCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc 000 18200 

C,..t.  00018300 

c«,.t.  oooieaoo 

C                     MAIN  LOOP  OF  PROGRAM  STARTS  HERE  00018500 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 000 18600 

C.,..,  00018700 

CMm.           CHECK   IF   TIME   TO  STORE  VALUES  FOR  PLOT  00018800 

Cimii  00018900 

250   IFtTIME.LT. TTPL)   GO  TO  35  00019000 

INDX  a  TIHE/PSTSZ  ♦   1  00019100 

T(INHX)    a   TIME  00019200 

Yi(l,INQX)   a  ASP  00019300 

Yi(2,INDX)   a  CON  00019400 

Y2(1#INPX)   a  ANN  00019500 

V2(2#IN0X)   a  PER  00019600 

Y3(1,INDX)    s   HERB  00019700 

Y3(2#INDX)   a  SHU  00019800 

Y4C1,INDX)   ■  ASP1  00019900 

Y4(2,INDX)   a  aSP2  00020000 

Y4(3,INDX)   a  A SP3  00020100 

Y5(1,INDX)   a  CON1  00020200 

Y5(2#IN0X)   a  C0N2  00020300 

Y5(3,INDX)   a  CON3  00020400 

CtfMf  00020500 

Cff9.tt           UPDATE  TTPL  FOR  NEXT  PLOT  TIME  00020600 

C«  •  #  • ,  00020700 

TTPL  a  TTPL   ♦  PSTSZ  00020800 

C.,,,.  00020900 

C                     CHECK   IF  TIME  TO  PRINT  00021000 

C...,,  00021100 

35  IF   (TIME.LT.TTPR)GO  TO  300  00021200 

C.tt,  00021300 

Ct....           WRITE  OUT  VARIOUS  VARIABLES  AND  TIME  00021400 

C   00021500 

WRITE  C6# 500) TIME  0  0021600 

500  FqRMATCIHOiIOX, 'TIME1 ,P7.2)  00021700 
WRITE (6,501 ) ASP. CON, SHU  0002180  0 

501  FORMAT ( 1  HO ,  '       ASP  a « , F20 , 6 , 5x # '  CON  a',F20,6,5x#  00021900 
/•       SHU  a»  ,F20.6)  00022000 

WRITE (6, 502) ANN, PER, HERB  00022100 

502  FORMAT ( 1 H   #'       ANN  a '  , F20 • 6 , 5x »  '  PER  a,,F20,6l5x#  00022200 
n     HPRB  a»F20,6)  00022300 

WRITE (6,507) ASPl. ASP?, A SP3  0  0  02240  0 

507  FORMAT  ( 1 H   ,«     ASPl   a^O^Sx.'  ASP2  a',F20.6,5X#  00022500 
/»     ASP3  a» ,F20.6)  00022600 

WRITE (6, 50 8) CON  1 ,CON?,CON3  00  02270  0 

508  FORMAT ( 1H   ,»     CQNJ   a • , F20 , 6 , 5x i 1  C0N2  e»,F20.6,5x*  00022800 
/■     CON3  ai,F20.6)  00022900 

WRITE (6,50 9) NOC ON, SUC #  SEED  0O023000 

509  FORMAT ( 1 H         NOCON  a'F20.6,5x,'  SUC  a»,F20,6,5X,  00023100 
/'     SEED  at,F20.b)  00023200 

WRITE (6,51 0)BI0SUM,SUCSUM#SEESUM  0  002330  0 

510  FORMAT ( 1H   ,'BIOSUM  a  •  F20 > , 5x t  '  SUC SUM  a',F20.6,5x,  00023400 
/'SEESUM  a«,F20.6)  00023500 

WRITE (6,51 1 ) AWPC  00023600 

511  F 0 R M A T ( 1 H  AWPC  a»,F?0.6)  O0O23700 
C,,..,  00023800 

C                       UPDATE   TTPR  FOR   NEXT   PRINT  TIME  00023900 

C....f  00024000 
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0  0  0  2  4  1  0  0 
00024200 

0002 a 3oo 

00024 U 00 
0002U500 
00024600 

BASIS,  0O024700 

00024600 
00^24900 
00025000 
0  0025100 
00025200 
0  0025300 
0OO25400 
00025500 
00025600 
00025700 

00025600 
00025900 
00026000 
00026100 
00026200 
00026300 
00026400 
00026500 
00026600 
00026700 
00026300 
00026900 

00027000 
0  0  0  2710  0 
00027200 
0  0  027300 
00027400 
00027500 
00027600 
00027700 
00027800 
00027900 


C.tt             COMPUTE  FLOW  1«         ASPI   GENERATION  00028000 

C,...,             AFFECTED  BY  ASP,   C0N#   SHU,    ASPI,   ASP2#   ASP3,   CONN  00028100 

C.....  00028200 

Zi«EXPC«P(33)*PASP)  0  0  028300 

Z2«EXP(*P(34)#PC0N)  0  00284  0  0 

Z3«l .♦CP(35)«1 , )*P3HU  0002850  0 

Z«»t .♦(P(36)-l t )*PASP1  0  002860  0 

Z5«l .♦(PC37)«1 ,)*PASP2  00028700 

Z6«1.*<PCJ8)»1.)*PASP3  0002880  0 

Z7»EXPC*PC39)*PC0NN)  0  0  02890  0 

FL0W14  »  P(40)*Zl*Z2*Z3*Z4*Z5*Z6*Z7  00029000 

Civiiff,  00029100 

Ct«ttt           COMPUTE  FLOW   15           ASPi   MORTALITY  00029200 

C,,.,           AFFECTED  BY  ASPI  00029300 

C,.,,f  00029400 

FL0W15  a  A8Pl*P(at)  00029500 

C,*tft  00029600 

C,!,.,            COMPUTE  FL0W16        ASPI  GRADUATION  00029700 

C,.,t,             AFFECTED  BY  ASP # CON # SHU #  ASPI t ASP2# A 5P3# CONN # HERB  00029800 

C,.,,.  00029900 

Zl«PU2)*(l.-P(42n*FxP(-PC43)*PASP)  O003ooo0 

Z2»EXP(«P(44)*PC0N)  0  0030  100 


TfPR  *  TTPR  ♦  STEP 

C  • 

C.!...  CHECK   IF   TIME  TO  STOP  RUN 

*300*IF(TEND,LEtTIME)   GO  TO  45 
C,f 

C..»tt  CALCULATE  PROPORTIONS  TO  PUT  ON  A  RELATIVE 

PSHUH  a  SHUH/MSHUH 
IFtPSHUH^GT. 1 . )PSHUH  s   t . 
PASPS  ■  ASPS/MASPS 
IF(PASPS.GT.i.)PASPSxi. 
PHERB  «  HERB/MHER8 
IF(PHERB.GT.l.)PHER8xl. 
PASP  ■  ASP/MASP 

IF(PASP.GT.1.)PASP»1. 
PCON  »  CON/MCON 

IF(PCONtGT.l.)PCON»l. 
PSHU  ■  SHU/M$HU 
IF(P$HU,GT.l.)PSHUsl. 
PASP1   ■  ASP1/MASP1 
lF(PASPl.GT.l.)PASPlsi. 
PASP2  a  ASP2/MASP2 
IF(PASP2.GT.l.)pASP2al. 
PASP3  «  ASP3/MASP3 
IF(PASP3.GT,l.)PASP3sl. 
PCONi  a  C0N1/MC0N1 
IFCPC0N1 .GT. 1 . )PC0N1»1 . 
PC0N2  a  C0N2/MC0N2 
IF(PC0N2.GT.i;)PC0N2al. 
PC0N3  a  C0N3/MC0N3 
IF(PC0N3.GT,1,)PC0N3»1. 
PCONN  a  CONN/MCONN 
IF  (PC0NN,CT.l.)PC0NN8l, 
PPER  a  PER/MPER 
IF(PPERtGT.l.)PPERal. 
PANN  a  ANN/MANN 
IF(PANN,GT.l.)PANNal. 
Ct.it.  CALCULATE  THE  VARIOUS  FLOwS 
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Z3»itt(P(45)«l.)*PSHU  0  0  030200 

Z«»i,f(PU6)»lf  )*PASP2  00030300 

ZSal.^(P(47)«l«)*PA8P3**P(4A)  00030400 

Z6sEXP(»P(49)*PCQNN)  Q0o3o500 

Z7«l .♦CPC50)»1 . )*PHERB  00030600 

FIQW16  a  ASPUP(5l)*Zl*Z2*Z3*Z4*Z5*Z6*Z7  00030700 

Ct9Vf  00030800 

C«.ft  COMPUTE   FL0W17         ASP2   MORTALITY  00030900 

Ci.mi  AFFECTED  BY  ASP2  00031000 

C,,m<  00031  100 

FL0W17  a  ASP2*P(52)  00031200 

C99.99.  00031300 

CMMi  COMPUTE  FL0W18         ASP2  GRADUATION  00031400 

Ct(iii  AFFECTED  BY  ASP,   CON,   SHU,    ASP2,   ASP3,   CONN  00031500 

Cv,»t«  00031600 

Zl«P(53)*Cl,»P(53))*ExP(-P(54)*PASP)  00031700 

Z2«EXPCPC55)*PC0N)  0  00  31800 

Z3«l.t(P(56)»t«)*PSHU  0  0031900 

Z4«l ,t(P(57)»l . )*PASP3  0  00  320  0  0 

Z5»EXP(«P(58) *PC0NN)  00032100 

FLOwlS  a  ASP2*P(59)*Zl*Z2*Z3*Z<i*Z5  00O32200 

C«««r«  00032300 

C«ffi«  COMPUTE  FL0W19  ASP3  MORTALITY  00032400 

C,....  AFFECTED  BY  ASP3  00032500 

C   00032600 

FL0W19  a  ASP3*P(60)  00032700 

C,99t,  00032800 

C9(.f,  COMPUTE  FLOW20  ASP3  GRADUATION  00032900 

C, . . t .  AFFECTED  BY  ASP #   CON,   SHU ,    A SP3  00033000 

C9.999  00033100 

Zl«P(61)f(l .•P(61) )*EXP(»P(62)*PASP)  0  00  3320  0 

Z2apXP(«P(63)*PC0N)  0  003330  0 

Z3al .♦(P(64).l . )*PSHU  0  0033400 

FLOW20  a  ASP3*PC65)*Z1*Z2*Z3  00033500 

Ctc999  00033600 

C   COMPUTE  FLOW  1         ASP3  GRADUATION  CONVERSION  00033700 

C.99.9  (ASP  GENERATION)  00033800 

C,,99,  00033900 

FLOW  1   a  FL0W20*P(66)  00034000 

C9t9l9  00034100 

C,,m.  COMPUTE  FLOW?  ASP  GROWTH  0003U200 

Cf9,tf  AFFECTED  BY  ASP,   CON  00034300 

C9.999  00034400 

Zl«l  ••(U«P (67) *,99/P (68)  )*PASP  0  003450  0 

Z2«lt*CP(69)»lt)*PCON**P(70)  00034600 

FL0W2  a  ASP*P(68)*ZUZ2  00034700 

C99ll9  00034800 

C9|999  COMPUTE   FL0W3  ASP   MORTALITY  00034900 

C,|9,9  AFFECTED  BY  ASP  00035000 

C,.999  00035100 

FL0W3  a  ASP*P(67)  00035200 

C99999  00035300 

C9|9t9  COMMRUTE   FL0W21  C0N1    REGENERATION  00035400 

C,,mi  AFFECTED  BY   ASP,   CON,   C0N2,   C0N3,   SHUH ,  HERB  00035500 

C99999  0O035600 

Zl»P(7i )♦( 1 9«P(7n )/P(7?)*PASP  000  35700 

IF   (Zi.GT.l.)   Zl=ll  00035800 

Z2al .♦(P(73)«l . )*PC0N**P(74)  000  35900 

Z3ai .♦(P(75)»l . )*PC0N2  0  00  360  0  0 

Z4ai , t(P (76) »1 , )*PC ON 3  00036100 

Z5«l .♦(P(77)»l . )*P*U"U  00  0  3620  0 
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Z6*i (P(78)«l . )*PHF_RB  00036300 

FL0w2l   s  P(79)*Z1  *Z?*Z3*Z  Wa*Z6  00036400 

C  •  .  f  t  •  0^036500 

C                      P(80)    YEAR   T I  Hp   DFLAV    In   CONFER  REPPODUCTION  00036600 

C  ,  .  ,  ,  ,  00036700 

IF    (TIME.LE.P(80))    Fl0*2ia0,  00036800 

C,,,,,  00036900 

C                         COMPUTE   FLOw?2           C0N1    MORTALITY  00037000 

C,,,M              AFFECTED   9Y   CON1  00037100 

Cf....  00037200 

Fl0w22  a  C0Nl*P(8n  00037300 

C  ,  ,  ,  ,  ,  0  0  03  7  400 

C,f.i.              COMPUTE   FL0W23           C0N1    GRADUATION  00037500 

C                         AFFECTED   BY   ASP,    C ON , C U N 1 , C r>N2 , C 0 N3 , SHUH , HERB  00037600 

C..,#.  00037700 

Zl*l .♦(P(82)»l .)*PA3P  0O037800 

Z2»l ,t(P(83)»l,)*PCQN  00037900 

Z3»l.*(P(84).l.)*PC0N2**P(85)  00038000 

Z^*l •♦<P(86)»1 ,)*PC0N3  00038100 

Z5»l .t(P(87)«l.)*PsHUH  00036200 

Z6»l  •♦(PC88)t»l  .  )*PHER8  00038300 

FLOW23  a  CON) ★PC89)#2l*Z2*Z3*Za*Z^*Z^  00038400 

C.m  0  00  3850  0 

C,f.»,            COMPUTE   FL0W24           cON2   MORTALITY  00038600 

C.,,..            AFFECTED   BY   CON?  00038700 

C..,.,  00038800 

FL0w24   a  C0N2*P(90)  00038900 

C,|M.  00039000 

c                 Compute  flowps     con2  graduation  00039100 

C,itt.             AFFECTED  9Y   ASP,C0N#C0N2,cnN3,  00039200 

C,.»..  00039300 

Zi»l .♦(P(91)»l . )*PASP  000  39400 

Z2*l , ♦ CP  C92 ) • 1 , )*PCON  00039500 

Z3sl,*(P(93)-l.)*PC0N2**PC94)  00  0  3960  0 

Za'l ■♦(p(95)»l. )*PC0N3**P(96)  000 39700 

FL0^25  a  CON?*PC973*Zl*Z2*Z3*Z4  00039800 

Cttitf  00039900 

C....t              COMPUTE   FL0WP6           CC1N3   MORTALITY  00040000 

C.tt              AFFECTED   9Y   C0N3  00040  1  00 

C   00040200 

FL0w26  a  C0N3*P(98)  00040300 

C,,,,,  00040400 

C..t».             COMPUTE   FL0w?7           CON3  GRADUATION    (N0C0N  GENERATION)  00040500 

C....,              AFFECTED   BY   ASP,    CON,    C0N3  00040600 

C   00040700 

Zl»l .♦(PC99)»1 . )*PASP  00040800 

Z.2»l.fCP(l00)«l.)*PCoN  0  0  04  090  0 

FU0W27  a  CON3*P(lon*ZUZ2  00041000 

C  «  .  .  t  •  0  0  04M00 

C...t.            COMPUTE   FL0W4           C0N3  GRADUATION   CONVERSION  0.0041200 

C...(.                                                 (CON  GENERATION)  00041300 

C.,tt.  00041400 

FL0W4   a  FL0^27*P  ( 1 0?)  00041500 

C...f  00041600 

C.....           COMPUTE   FLOW5           CON  GROWTH  00041700 

C...            AFFECTED   BY   ASP,    CON  00041800 

C,  ,  ,  ,  ,  00041900 

Zlapt 103) + ( 1 .»P( 103) ) *EXP(.P( 104) *PASP)  00042000 

Z2*l 1 ,*P( 105)*,99/P( 106) )*PCHN  000  42100 

FL0W5  a  CDN*P(  lok^7,W3  0OQ42200 

C.  •  «  •  .  00042300 
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C,..f,             COMPUTE                          NOCON  MORTALITY  00042400 

C.?.«.           AFFECTED  BY     NOCON  00042500 

Cn.fi.  00042600 

FLOwSO  s  N0CQN*P(105)  00042700 

C.,,.  00042800 

C                       COMPUTE  FLOWfc           C ON  MORTALITY  0O0429O0 

C....,           AFFECTED  8Y  CON  00043000 

C...a,  00043100 

FL0W6  a  CON*P(105)  00043200 

C.fftt  00043300 

C,,tM              COMPUTE   FL0W7           ANN  PRODUCTION  00043400 

CtfiM           AFFECTED  BY  ASP ,   C0N#   PER »   SHUH.CONN  00043500 

C,...,  000^3600 

Zl   s  EXP(»P(107)*PASP)  00043700 

Z2  a  P(lOfl)f  Cl.«»PCl08n*F*PC-PMo9)*PCnKj)  00043800 

Z3  a  1,  +  CP(U0)«1, UPPER  00043900 

Z4  s  EXP(«P(U1)*PSHUH)  00044000 

Z5  a  P(ll2)t(l,»P(tl2))*ExP(-P(lt^)*PC0NN)  00044100 

FL0W7  8  P(114)*Z1*Z2*Z3*Z4*Z5  00044200 

C,,,t.  00044300 

Ctfttf            ANNUAL   MORTALITY   ASSUMED   TO  HE    100%  EACH   YEAR  00044400 

C«|M,  00044500 

FLOWS   »   ANN/DT  00044600 

C, 9ta,  00044700 

C»,.i,           COMPUTE  FL0W9           PF«  PRODUCTION  00044800 

Cifi!f           AFFECTED  BY  ASP,   CON,   ANN,   SHUH.   CONN  00044900 

C.,.«  00045000 

Zl   a   l.f(P(115)«l.)*PASP  00045100 

Z2  a  P(116)+(1.«P(U6) )*FXP(«P( I 17)*PC0N)  00045200 

Z3  8   l.+(P(118)»l,)*PANN  00045300 

Z4  a   i,*(PC119)wl.)*PSHUH  00045400 

Z5  a  P(120)*(U«P(120))*PXP(-P(121)*PCONN)  00045500 

FL0W9  s  P(122)*Z1#Z2*Z3*Z4*Z5  00045600 

00045700 

C..#tt  PER  MORTALITY   (A80vF  GROUND)    ASSUMED  TO  BE   100*  EACH  YEAR  00045800 

C   00045900 

FL0W10  s  PER/DT  00046000 

C99 »t>  00046100 

C9Mtf            COMPUTE   FL0W11            SHU  REGENERATION  00046200 

C                     AFFECTED  BY  ASP , CON • SHU  $    ASPS,   CONN  #  HERB  00046300 

C,.,,.  00046400 

ZlaE)(P(«P(  123)*PASP)  00046500 

Z2beXP(i.P(  124)  *PC0N)  000  4660  0 

Z3  a   l,t(PU25)-lt)*PSHU**P(1261  00046700 

Z4  s   l.f(P(127)*t.)*PASPS  00046800 

Z5  a  EXP(«P(128)*PC0NN)  00046900 

Z6  a   l.t(P(129)«i,)*PHERB  00047000 

FLOwli  s  P( i30)*Zl*z2*Z3*Z4*z5*Z^  00047100 

C.,.t.  00047200 

C9,tt«           COMPUTE  FLOW12           SHU  GROWTH  00047300 

Cf|tM           AFFECTED  BY  ASP,   CON*   SHU,    ASPS.   CONN,   HERB  00047400 

Ct .,t,  00047500 

Zl   a  1 ,*(P( 131 )-l , )*PASP  00047600 

Z2  s  1 ,*(P( 132)«t . )*PCOM**p( 133)  00047700 

Z3  s   1 t«(l ,»P(134)*.99/P( 135) )*PSHU  00047800 

Z4   8    1 ,f(P( 136)»1 . )*PASPS  00047900 

Z5  a   1 ,4(P( 137)-! . )*PCONN  00048000 

Zb  a   l.t(P(138)»t.)*PHERH  00048100 

FL0W12  a  SHU*P(  1  35 )  *  Z 1  *  Z?*  73*  Za  *  Z  W<>  00048200 

C   00048300 

C...«.            COMPUTE   FLOW13           SHU   Hf)R  T  A|  I  TV  00048400 
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C.,,,.  AFFECTED  BY  SHU  00048500 

C.....  000U86OO 

FLQwl3  =  SHU*P(l3a)  00048700 

C.,.,.  00048800 

C   UPDATE   STATE   VARIABLES  00048900 

C9.««.  00049000 

ASP  »  ASP  ♦   OT   *    (FLOWl    ♦  FL0W2  -  FLOWS)  00049100 

CON  s  CqN  ♦  OT  *    (FLO^   ♦  FL0*5  -  FL0^6)  00049200 

ANN   3   ANN   ♦   DT   *    (FL0W7   •   FLOWft)  00049300 

PEW   s  PEP   ♦   DT   *    (FL0W9  •  FL0W10)  00049400 

SHU  a  SHU  ♦  OT  *   CFLQWll   +  FL0W12  *  FL0*13)  00049500 

BIOSHR  s  BIQSQ»   f  OT   *   (.FL0Wl«rLOH2-FLOwa-FL0W5-  00049600 

lFL0w7itFl.Ow9.FLC1w1  1-FLOwi?)  00049700 

BIOSIN  n  BIOSIN  ♦  DT  *   ( F  LcH  3 1 F  L0W6  t  Flo*  9  *F  L  0*  1  0  ♦FLO*  1  3 )  00049800 

ASPl    s   ASPl    ♦   DT   *    (FL0W14   «  FLOW  15  •  FL0>'16)  0O049900 

ASP?  s   ASP?   ♦  OT   *    (FLOwtb  -   FL0*l7   m  FL0W13)  00050000 

ASP3  »  ASP3   t  DT   *   C  FLOW  1 8  m  FL0*i9  -  Flq«2o)  00050100 

SUCSOR  s  SUCSOR  +  OT   *   (•FL0W14)  00050200 

SUCSIN  s  SUCSIN  +  DT  *   ('FLOW  1 5  +  FLOw  1 7  +  FLCW 19+FLOW20 )  00050300 

CONl   s  C0N1   ♦   Dj   *    (FLU^2l   -  FLDW?2  •  FL0^23)  00050400 

C0N2  *  CON 2  ♦  DT  *   (FL0W23  -  FLOw?a  •  FL0W25)  00050500 

C0N3  »  CONS  ♦  DT   *    (FL0w?5  •  FLOw?fe  •  FL0w27)  00050600 

SEESOR  a  SEE30P  ♦  OJ  *   C»FL0*2i)  00050700 

SEESIN  c  SEESIN  ♦  DT  *   (FLOW22+FLOW24+FL0 w26*FL0W27 )  00050800 

CtVVf.  00050900 

C,..i.  UPDATE   VARIOUS   AUXILIARY   VARIABLES  00051000 

Cmim  00051  100 

HERB  »  ANN  ♦  PER  00051200 

AS^S  B  P(27)*ASP1   ♦  P(?8)*ASP2  ♦  P(?9)*ASP3  00051300 

SHUH  a  SHU  +   ASPS  00051400 

CONN  a  p(30)*CONl    ♦  P(3l)*C0N2  ♦  P(3?)*CON3  00051500 

NOCoN  a  NQCqn  ♦  Ot  *   CFLO*27  -FlqwSO)  00051600 

SUC  a  ASPl    +  ASP2  *  ASPS  00051700 

SEED  a  CONl   f  C0N2  ♦  CONS  00051800 

8I0SUM  »  BIQSOR  ♦  ASP  ♦  CON  t  a-jn  ♦  PFR  +  SHU  +  BIOSIN  00051900 

SUCSUM  9  SUCSOR  ♦  ASPl   +  A3P2  ♦  ASP3  +  SUCSIN  0OO52000 

SEESUH  3  SEESOR  ♦  CONl   +  C0N2  ♦  CHN3  *  SEESIN  0O052100 

IFCnqCqn.GT.0.0)   Go  TO  240  00052200 

AWPCa.l  ,  0O0523O0 

GO  TO  242  00052400 

240  AWPCaCON/NOCON  00052500 

242  continue  ooo526oo 

C   00052700 

C   UPDATE  T I  HE  00052800 

C,,,,,  00052900 

2-49     TIME     a  TIME     ♦  DT  00053000 

C   00053100 

GO  TO  250  00053200 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCrCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc 0  0 05330  0 

C...,,  MAIN  LOOP  OF  PROGRAM  ENDS  HERE  00053400 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 00053500 

C.,,t.  00053600 

C   DR Art  GRAPHS  OF  SELECTED  STATE  VARIABLES  00053700 

C»....  00053800 

45  WRITE(6,512)  00053900 

512  F0RMATClHl,30Xf  •  1  o  &5P  2  *  Con',/)  00054000 
CALL  PL0TSC2,  INDX#T,  YD  00054100 
rtRlTE (6,513)  0OO54200 

513  FQRMAT( 1 H 1 ,30X, » 1  8  ANN  2  a  Pprpi,/)  00054300 
CALL  PL0TS(2#INDX,T,Y2)  00054400 
WRITE(6,514)  0005a5o0 
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APPENDIX  A  (cont.) 


5ia  FORMAT(1H1,30X, • 1   a  HERB           2  a  SHU',/)  00054601 

CALL  PL0TSC2, INDX, T, Y3)  oo<)54700 

WRITE (6,515)  0005U800 

515  FORMAT  ( 1H1 ,  30X»  '  1  a  ASPl  2  s  ASP?  3  a  ASP  3 '  ,  / )  0005/4000 
CALL  PL0TS(3# INDX,T, Y4)  00055000 
HRITE(6,516)  00055100 

516  FORMAT(lHl,30X, • 1  a  CQNl  2  a  C0N2  3  a  C0N3' ,/)  00055200 
CALL  PL0TS(3# IMDX#T,Y5)  00055300 

END  00055400 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc 00 05550  0 

C...           END  OF  MAIN  PROGRAM  00055600 

ccccccccccccccccccccccccccccccccccccrcrrccr:ccccccccccccccccccccrcccccccrooo557oo 

C   00055800 

C,,.tt  00055900 

Ci.mi           GRAPHING  SUBROUTINES  00056000 

cccccccccccccccccccccccccccccccccccccrcrcccccccccccccccccccccrcccccccccrooo56ioo 

SUBROUTINE  PLOTS ( NPLOTS ,  NPT S # X # Y )  00056200 

C0MM0N/I0/NI ,N0  00  05630  0 
DIMENSION  X(NpTS)#Y(NPLOT$,NPTS)#xSCl  ( 1 ®  I)  t YSCL ( 52 ) , I L I NE ( 1 0 ! )  00056400 

NI   3  5  00056500 

NO  a  6  00056600 

XMlNaX(l)  00056700 

XMAXaXC 1 )  00056800 

DO    1    IaJ,NPTS  00056900 

lF(X(I),GT.XMAX)XMAXaX(I)  00  057000 

1  IF(X(I)tLT.XMIN)XMlNaXm  00  057100 
YMlNaY(Ul)  00057200 
YMAX»Y(1#1)  00057300 
00  2  InWNPLOTS  00057400 
DO  2  J"1#NPT$  00057500 
IFCYdf  J)  ,LTf  YMIN)YHlNaYf  It  J)  00  057600 

2  IF (Y (I# J) ,GT. YMAX) YHAXaY (1 1  J)  0  0  05770  0 
7           SCALX«(XMAX«XMIN)/100.  00057800 

SCALy«(YMaX«YMIn)/50I  0  0  05790  0 

00  3   I»l»101  00058000 

XSCL(I)axMINfFLOAT(I»l)*SCALX  0  0  058100 

IF(I.GT.52)G0  TO  3  00058200 

YSCL (I ) a YMAX. FLO AT (1.1 )*SCALV  0  0  05830  0 

3  CONTINUE  00058400 
C..,..  00058500 

CALL  GZROCO, ILINE)  00058600 

C.,ft,  00058700 

DO  5  Kal,51  00058800 

call  gzro(k# iline)  00058900 

call  skpt(nplots#npts#k. ** yscl,xmINiscalx, iline)  00059000 

5  call  typit(k# yscl, iline)  00059100 
call  g2r0(0# iline)  00059200 

WRITE(N0#6) (X$CL( I) , ial , J01 ,20)  0*059300 

6  FURMAT(17X,F7,2,5(13X,F7.2) )  0  0  0594  0  0 
RETURN  00059500 
END  00059600 

C...  00059700 

SUBROUTINE   SKPT(NPLOTS#NpTS,K#X#Y, YSCL*  XM IN,  SCALX# ILINE)  n 0  05980  0 

C,....  00059900 

DIMENSION  X(NPTS) #Y (NPLOTS, NPTS) »  TL INE ( 1 ) # YSCL ( 1 ) #  ISYM (  1  1 )  00060000 

DATA   ISYM/1H1 , 1 H2, 1H3, 1H4, 1H5, 1 M6# 1H7, IH8, 1H9, 1  HO, 1HC/  0  0  060  100 

DATA   IBLK/1H  /  00060200 

DO   1    I  a  1 , NPLOTS  00060300 

DO   1   J«1#NPT$  00060400 

IF(Y(I#J),GT.YSCL(K))GO  TO  1  OOO6050O 

IF(Y(I#J),LE.YSCL(Kfl))GO  TO   1  00060600 
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APPENDIX  A  (cont.) 


Mx I  NT (  (X(J)*XMIN)/SCALa,M  tbj 

u  u  u  0  v  f  n  u 

Ir  (M ,UT , 1 J 60   TO  1 

U  v  0  Q  U  O  0  U 

IF  (M.GT, 101  )G0  TO  1 

0  U  (1  QU  90  u 

IF ( ILINE  (M)  ,EQ,  IBLK  )G0  TO  ? 

00 u ©1000 

ILINc  In)  *ISY"1  (  1 1  J 

u  n  u  0  1  1  0  u 

GOTO  1 

A  A  A  Lm.  A    T  *  A 

OOuOjcOO 

2 

I L I N  E (M)bISYH ( I ) 

A    A    ii           4     7    4  A 

OOOol  iO'' 

1 

n   ^%  »  |           Skill  4k- 

CONTINUE 

A  f\  r\         4    /i  4  A 

RETURN 

0  0  0  0  1  5  0  0 

END 

00  0O1700 

SUBROUTINE  TYPIT(K#Y5CL#ILINtJ 

AAAK1 ftnA 

m 

A  A  A  <«  1  Q  A  A 

O'lOOj  "1)0 

CQnMQN/IO/NIfNO 

AAA*h3AaA 

DIMENSION   Y5C.L  1 1  J  #  ILINc  (  1  J 

A  A  A  ^  3  1  A  A 

f  m  mm  IS  a.  il 

JS'K  ▼  <4 

aAAA«33aA 

0  u uoecO  0 

I r C ( J5» JS/5*5 ) t NE , 0 ) GO  TO  1 

1 J  u  1 1  ft  ^ JUU 

WKITE(N0#6)  YSCL(K) 

AAA*»3/JAA 

6 

FORMAT C 1  0x# E 1 0 t 3# 2H  ♦) 

i.i  mm.   m>    _  V*    A    .   .    _           ■%   \     /    —   ,       «■       ,  m     *           a                        •           4     A    4  \ 

WRITECNQ#2) ( I  LINE  Cl)#I«l#*0lJ 

a  a  n  t  T  L  *\A 

2 

F0RMAT(1HT#20X# 101Ai»2X#?H  +  n 

AAA  /  "S»»AA 

0  0  u  0  c  7  0  " 

GO  TO  3 

U  ' ' '  j  O  c  O  0  U 

• 
1 

WRITE (N0#5) 

AAAAi3QAA 

5 

format (2ix0 ih«j 

A  n  m  L  7  A  a  n 

0  0  ()  O  3  0  0  U 

WRlTE(N0#4)(lLlNE(I),Ial,10n 

A  A  n      1  1  aA 

000031  no 

FORMAT(lHf#20x#10lAlf2X#2H*I) 

A  A  A  ^  1  ^  A  A 

3 

RETURN 

0  0063300 

kl  l\ 

END 

/I  ft  Al  7  /t  a  A 

1)0063*100 

00  0&3500 

SUBROUTINE  GZRO ( K # I L I NE ) 

AAA*»1i»AA 

0  u  0  O  3  0  0  0 

f 

0  0  0  0  3  f  0  U 

COMMON/IO/NI #  NO 

aaa^IAaA 

DIMENSION  ILINE(I) 

AAA*k"lOAA 

0  0  0  6  3^00 

DATA   IPER# IN,  IBLK/1H. - 1H1 9 1H  / 

A  A  A  <•  /I  A  A  A 

(1  n  0  O  (4  0  0  0 

IF(K.E0,0)G0  TO  2 

A  A  A      /I  1  A  A 

0  O  (J  0  4  1  0  U 

DO  1   I»l , 101 

aaaa>/i3aA 
0  U  0  O  4  c  0  0 

i 

l 

IL I NE ( I ) si BLK 

AAAA/l^AA 

1)  0  1]  O  ■*  j  0  U 

RETURN 

A  A  A  it  /I  /I  A  A 

2 

DO  3  1*1*101 

AAAAtl^AA 

y  u  u  0  «*  j  (j  if 

3 

ILINE ( I  )»IPER 

aaaku4»aa 

P\                        It               •     SB     4                 J       A      A                 4  A 

Do  «  IS1#101#10 

A  A  l\  4*  /I  7  A  A 

ILINE(I)»IN 

AAAAZiAAA 
II  U  U  O  M  O  0  U 

WRITE(NQ#5) (ILINE(I) # 1*1  - 101 ) 

aaa*«/iOaA 

0 ')  0  0  n  y  0  'J 

5 

F0RMATCl9X,2H,.,l0lAl#2H.,) 

A  A  A  ^»  C  A  A  A 

RETURN 

000651 00 

END 

00065200 
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APPENDIX  A  (cont.) 
1  DATA 

i,,400,,lOO,,4, 

200000,, 250  0  0  0., 1000  0,, 3500. ,500, 

1  0  0  0  0  0  ,  ,  5  0  0  0  0  ,  ,  1 5  0  0  0  .  ,  1  0  0  0  0  0  .  t  1  0  0  0  0  , 

2500, , 1500  0 , # 1500  0. , 150  00, . . 1 

,l,,l,.l#500tf .1 

vl,«l,«l#»l#*l 

0,  , ,00b,  ,1 # ,8,  ,02 

,  a  ,  i  ,  a  ,  5  •  ,  6  ,  ,  .  8 

,3,, i#, 001, 6, ,10 oooo. 

,45,, 1,5, ,5,, ,5 
,  0  0  1 ,  .  0  0 1  ,  .  6  7  #  S ,  •  ,  6 
,25, ,25, , 1,5, ,5, 
,7, ,001, 5., ,2, ,1 
,1,5,,5,,,$,,15 
a,,  ,03, ,08, ,001, 1.2 
, 6,, 25, ,1,3, ,.2 
,  2 ,  ,  2 ,  ,  0  0  l ,  1  0  0  0  , ,  5 . 
,  7 , ,5, ,5, ,5,3, 
,2# ,7, ,7, .3, ,45 
,5,,3,,7,3t,i5 
3, « ,25, , 1 , ,5, , 3 
,2,6,,,5,3,,«02 
.07,2,  ,  ,005,6,  ,  .to 
2,  ,  ,005,6, ,500 . , ,5 
,007,&,# ,7, ,5, .007 
6, •  3500  , , a, , 6,  • , 1 
3,5, ,^,6,# ,8, 10  0. 
,5, .001 , ,5, ,03, ,2 
,5, ,01, ,4,0,0 
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APPENDIX  C 


Parameter  Names,  Definitions,  Values,  Range  of 
Values,  and  Units  for  the  ASPEN  Program. 


Param- 
eter 
number 

FORTRAN 
name 

Definition 

Value 

Range 

Units 

Source 

1 

MASP 

Maximum  value  for  aspen 
hiomass 

200,000 

kg/ha 

Zimmermann  1979  Rodin 

and  Ra7ilp\/irh  1QA7 

2 

MCON 

Maximum  value  for  conifer 

U\\J\  1  laoo 

250,000 

kg/ha 

Zimmermann  1979 

3 

MSHU 

Maximum  value  for  shrub 
biomass 

10,000 

kg/ha 

Unpublished  data 
(Bartos,  files) 

4 

MPER 

Maximum  value  for  peren- 
nial biomass 

3,500 

kg/ha 

Youngblood  and 
Mueggler  1981 

5 

MANN 

Maximum  value  for  annual 
biomass 

500 

kg/ha 

Youngblood  and 
Mueggler  1981 

6 

MASP1 

Maximum  value  for  ASP1 
reproduction  class 

100,000 

80,000-120,000 

Numbers/ha 

Baker  1925  and 
unpublished  data  (Bartos, 
files) 

7 

MASP2 

Maximum  value  for  ASP2 
reproduction  class 

50,000 

40,000-60,000 

Numbers/ha 

Unpublished  data 
(Bartos,  files) 

8 

MASP3 

Maximum  value  for  ASP3 
reproduction  class 

15,000 

10,000-20,000 

Numbers/ha 

Mueggler  and  Bartos 
1977 

9 

MCON1 

Maximum  value  for  CON1 
reproduction  class 

100,000 

80,000-120,000 

Numbers/ha 

Noble  and  Ronco  1978 

10 

MCON2 

Maximum  value  for  CON2 
reproduction  class 

10,000 

5,000-15,000 

Numbers/ha 

INT-1751  estimate 

i  i 

1 VI  uuINO 

MaYimnm  v/alit^  for  f^ON^ 

IVI dA  1 1  IIU1II   VdlUC   1  \Ji    UVIN J 

reproduction  class 

9  son 

1  noo-5  oof) 

Ni  imhpr«i/ha 

iiui  i  iuci  o'l  ia 

(Ml      I  1  \J  \   Co  1 1 1  1  Id  LC 

12 

MASPS 

Maximum  biomass  value 
for  the  sum  of  all  asDen 

1        1       III  w     \J  W  1  1  1     >■/  1     WAI  1     (AW       W  I  I 

suckers 

15,000 

10,000-20,000 

kg/ha 

INT-1751  estimate 

13 

1  \J 

MSHUH 

Mavimum  hioma^^  value 

IVIUAM  1  IUM  I    UIUI  1  luuu  VUIUU 

for  the  sum  of  all  aspen 
suckers  and  shrub  biomass 

15  000 

10  000-20  000 

kg/ha 

INT-1751  estimate 

II1!    1         liwl     C  JIM  1  1  Q  I  V> 

14 

MCONN 

Maximum  biomass  value 
for  the  sum  of  all  conifer 

eppH  1  i  nnc 

OCCU  11  1  1  y  O 

15,000 

10,000-20,000 

kg/ha 

INT-1751  estimate 

15 

ASP 

Initial  value  for  state 
variable  aspen 

0.1 

0-1. 

kg/ha 

—  — 

16 

CON 

Initial  value  for  state 
variable  conifer 

.1 

0-1. 

kg/ha 

17 

SHU 

Initial  value  for  state 
variable  shrubs 

.1 

0-1. 

kg/ha 

18 

PER 

Initial  value  for  state 
variable  perennials 

.1 

0-1. 

kg/ha 
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Param- 


eter 
number 

FORTRAN 
run  i  r\M vi 

name 

Definition 

Value 

Range 

Units 

Source 

19 

ANN 

Initial  value  for  state 
variable  annuals 

500 

250-500 

kg/ha 

Bartos  and  Mueggler 
1981 

20 

ASP1 

Initial  value  for  state 
variable  ASP1 

0.1 

0-1. 

Numbers/ha 

OH 

A  C  DO 

Abrii 

I  r*»  i  +  i  o  I    \  /  o  1 1  m   t/"\  i*   plot  e\ 

initial  vaiue  ior  siaie 
variable  ASP2 

.  i 

iNumDers/na 

22 

ASP3 

Initial  value  for  state 
variable  ASP3 

.1 

0-1. 

Numbers/ha 

23 

CON1 

Initial  value  for  state 
variable  C0N1 

.1 

0-1. 

Numbers/ha 

24 

CON2 

Initial  value  for  state 
variable  CON2 

.1 

0-1. 

Numbers/ha 

25 

CON3 

Initial  value  for  state 
variable  CON3 

.1 

0-1. 

Numbers/ha 

26 

NOCON 

Initial  value  for  conifer 
trees  on  site 

0 

Trees/ha 

27 

Factor  to  convert  sucker 
numbers  (ASP1)  to  biomass 

.006 

.004-. 006 

kg/sucker 

Bartos  and  Johnston 
1978 

28 

Factor  to  convert  sucker 
numbers  (ASP2)  to  biomass 

.1 

.04-. 11 

kg/sucker 

Bartos  and  Johnston 
1978 

29 

Factor  to  convert  sucker 
numbers  (ASP3)  to  biomass 

.8 

.6-1. 

kg/sucker 

Bartos  and  Johnston 
1978 

30 

Factor  to  convert  seedling 
numbers  (CON1)  to 
biomass 

.02 

.01 -.03 

kg/seedling 

Unpublished  data 
(Bartos,  files) 

31 

Factor  to  convert  seedling 
numbers  (CON2)  to 
biomass 

.4 

.2-.6 

kg/seedling 

Unpublished  data 
(Bartos,  files) 

32 

Factor  to  convert  seedling 
numbers  (CON3)  to 
biomass 

2.4 

2.-3. 

kg/seedling 

Unpublished  data 
(Bartos,  files) 

33 

Aspen  restriction  on  ASP1 
regeneration 

5. 

3.-7. 

Dimension- 
less 

INT-1751  Estimate 

34 

Conifer  restriction  on  ASP1 
regeneration 

6. 

4.-8. 

Dimension- 
less 

INT-1751  Estimate 

35 

Shrub  restriction  on  ASP1 
regeneration 

0.8 

.6-.9 

Dimension- 
less 

INT-1751  Estimate 

36 

ASP1  restriction  on  ASP1 
regeneration 

.3 

.1-.5 

Dimension- 
less 

INT-1751  Estimate 

37 

ASP2  restriction  on  ASP1 
regeneration 

.1 

0-.3 

Dimension- 
less 

INT-1751  Estimate 
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Param- 
eter FORTRAN 


number 

name 

Definition 

Value 

Range 

Units 

Source 

38 

ASP3  restriction  on  ASP1 
regeneration 

.001 

0-.2 

Dimension- 
less 

INT-1751  Estimate 

39 

CONN  restriction  on  ASP1 
regeneration 

6. 

4.-8. 

Dimension- 
less 

INT-1751  estimate 

40 

Unrestricted  ASP1  regener- 
ation rate 

100,000 

80,000-120,000 

Suckers/ha/ 

yr 

Baker  1925  and 
unpublished  data  (Bartos, 
files) 

41 

ASP1  mortality 

0.45 

.3-.6 

Yri 

Unpublished  data  (Bartos, 
files) 

42 

Aspen  restriction  on  ASP1 
graduation 

.1 

0-.2 

Dimension- 
less 

INT-1751  estimate 

43 

Aspen  restriction  on  ASP1 
graduation 

5. 

3.-7. 

Dimension- 
less 

INT-1751  estimate 

44 

Conifer  restriction  on  ASP1 
graduation 

5. 

3.-7. 

Dimension- 
less 

INT-1751  estimate 

45 

Shrub  restriction  on  ASP1 
graduation 

.5 

.3-7 

Dimension- 
less 

INT-1751  estimate 

46 

ASP2  restriction  on  ASP1 
graduation 

.001 

0-.2 

Dimension- 
less 

INT-1751  estimate 

47 

ASP3  restriction  on  ASP1 
graduation 

.001 

0-.2 

Dimension- 
less 

INT-1751  estimate 

48 

ASP3  restriction  on  ASP1 
graduation 

.67 

.33-1. 

Dimension- 
less 

INT-1751  estimate 

49 

CONN  restriction  on  ASP1 
graduation 

5. 

3.-7. 

Dimension- 
less 

INT-1751  estimate 

50 

Herb  restriction  on  ASP1 
graduation 

.8 

.7-.9 

Dimension- 
less 

INT-1751  estimate 

51 

Unrestricted  ASP1  gradua- 
tion rate 

.25 

.15-.35 

Yr-1 

INT-1751  estimate 

52 

ASP2  mortality 

.25 

.15-.35 

Yr~1 

Unpublished  data  (Bartos, 
files) 

53 

Aspen  restriction  on  ASP2 
graduation 

.1 

0-.2 

Dimension- 
less 

INT-1751  estimate 

54 

Aspen  restriction  on  ASP2 
graduation 

5. 

3.-7. 

Dimension- 
less 

INT-1751  estimate 

55 

Conifer  restriction  on  ASP2 
graduation 

5. 

3.-7. 

Dimension- 
less 

INT-1751  estimate 

56 

Shrub  restriction  on  ASP2 
graduation 

.7 

.5-.9 

Dimension- 
less 

INT-1751  estimate 

57 

ASP3  restriction  on  ASP2 
graduation 

.001 

0-.2 

Dimension- 
less 

INT-1751  estimate 
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Param- 
eter FORTRAN 


number 

name 

Definition 

Value 

Range 

Units 

Source 

58 

CONN  restriction  on  ASP2 

nraH  i  lat  i  rin 
y  1  aU  Uau  Ui  l 

5. 

3.-7. 

Dimension- 

coo 

INT-1751  estimate 

59 

Unrestricted  ASP2  gradua- 
tion r?»tp 

1 1       1  1  O  IC 

.2 

.1-.3 

Yr~1 

INT-1751  estimate 

60 

ASP3  mortality 

.1 

.05-. 15 

Yr"1 

Unpublished  data  (Bartos, 
files) 

61 

Aspen  restriction  on  ASP3 
nradi  latinn 

.1 

0-.2 

Dimension- 
|p^^ 

ICOO 

INT-1751  estimate 

62 

Aspen  restriction  on  ASP3 
nraduation 

5. 

3.-7. 

Dimension- 
less 

INT-1751  estimate 

63 

Conifer  restriction  on  ASP3 

n raH i i  a  t  i  n  n 

5. 

3.-7. 

Dimension- 

I  t?oo 

INT-1751  estimate 

64 

Shrub  restriction  on  ASP3 

n raH i  la  t  i  n  n 
y i  au uduui  i 

.8 

.6-.9 

Dimension- 

Incc 
i  coo 

INT-1751  estimate 

1  InroctriptoH  AQP^  nraHna- 

tion  rate 

1*5 

1  -  7 

Yr"1 

I  MT-1  7  m  oct  i  ma  to 

ACD'i  nrarli ifltinn  rnnvpr- 

no  ui  uuuouui  i  uui  ivci 

sion  to  biomass 

A 

3  -5 

k  n/t  rpp 

r\ y /  Li  cc 

Rartnc  anri  IrihnQtnn 
uai  iuo  diiu  JwiiiiolUll 

1978 

67 

A<?npn  mortality 

.03 

.02-.04 

Yr"1 

INT-1751  pctimatp 

68 

Unrestricted  aspen  growth 
rate 

.08 

.06-1. 

Yr~1 

INT-1751  estimate 

69 

Conifer  restriction  on 

00      1 1  y  i  <j  w  l  i  i 

.001 

0-.02 

Dimension- 

INT-1751  estimate 

70 

Conifer  restriction  on 
a^npn  nrowth 

UOp&l  1    UI  Ull  LI  1 

1.2 

1.-1.4 

Dimension- 
less 

INT-1751  estimate 

71 

Aspen  restriction  on  C0N1 

rononcrat  inn 

.6 

.4-.8 

Dimension- 

ICOO 

INT-1751  estimate 

72 

Aspen  restriction  on  C0N1 

rpnpnprfl  tion 
1  CUgl  ICI  allUI  1 

.25 

.15-.35 

Dimension- 

ICOO 

INT-1751  estimate 

73 

Conifer  restriction  on 

v-'^/iN  i  i cy  t?i  ici  duui  i 

.1 

0-.2 

Dimension- 

I  PC  c 
Icoo 

INT-1751  estimate 

74 

Conifer  restriction  on 

HON1  rpnpnpration 

^  i  \  '  1  1  1     1  Cy  CI  lu  1  O  LI  Wl  1 

3. 

1.5-5. 

Dimension- 
less 

INT-1751  estimate 

75 

C0N2  restriction  on  C0N1 
regeneration 

.2 

0-.4 

Dimension- 
less 

INT-1751  estimate 

76 

CON3  restriction  on  CON1 
regeneration 

.2 

0-.4 

Dimension- 
less 

INT-1751  estimate 

77 

Shrub  and  aspen  sucker 
restriction  on  C0N1 
regeneration 

.2 

0-.4 

Dimension- 
less 

INT-1751  estimate 

78 

Perennial  and  annual 
restriction  on  C0N1 
regeneration 

.001 

0-.2 

Dimension- 
less 

INT-1751  estimate 
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Param- 
eter FORTRAN 


number 

name 

Definition 

Value 

Range 

Units 

Source 

79 

Unrestricted  CON1 
regeneration  rate 

1,000 

750-1,250 

Seedlings/ 
ha/yr 

Noble  and  Ronco  1978 

80 

Time  delay  in  years  for  con- 
ifer reproduction 

5. 

3.-10. 

Yr 

INT-1751  estimate 

81 

CON1  mortality 

.7 

.5-.9 

Yr-1 

Noble  and  Ronco  1978 

82 

Aspen  restriction  on  CON1 
graduation 

.5 

.3-.7 

Dimension- 
less 

INT-1751  estimate 

83 

Conifer  restriction  on 
CON1  graduation 

.3 

.1-.5 

Dimension- 
less 

INT-1751  estimate 

84 

CON2  restriction  on  CON1 
graduation 

.5 

.3-7 

Dimension- 
less 

INT-1751  estimate 

85 

CON2  restriction  on  CON1 
graduation 

3. 

1.-4. 

Dimension- 
less 

INT-1751  estimate 

86 

CON3  restriction  on  CON1 
graduation 

.2 

0-.4 

Dimension- 
less 

INT-1751  estimate 

87 

Shrub  and  aspen  sucker 
restriction  on  CON1 
graduation 

.7 

.5-.9 

Dimension- 
less 

INT-1751  estimate 

88 

Perennial  and  annual 
restriction  on  CON1 
graduation 

.7 

.5-.9 

Dimension- 
less 

INT-1751  estimate 

89 

Unrestricted  CON1 
graduation  rate 

.3 

.1-.4 

Yr"1 

INT-1751  estimate 

90 

CON2  mortality 

.45 

.3-.6 

Yr"1 

Noble  and  Ronco  1978 

91 

Aspen  restriction  on  CON2 
graduation 

.5 

.3-7 

Dimension- 
less 

INT-1751  estimate 

92 

Conifer  restriction  on 
CON2  graduation 

.3 

.1-.5 

Dimension- 
less 

INT-1751  estimate 

93 

CON2  restriction  on  CON2 
graduation 

.7 

,5-.9 

Dimension- 
less 

INT-1751  estimate 

94 

CON2  restriction  on  CON2 
graduation 

3. 

1.-4. 

Dimension- 
less 

INT-1751  estimate 

95 

CON3  restriction  on  CON2 
graduation 

.5 

.3-7 

Dimension- 
less 

INT-1751  estimate 

96 

CON3  restriction  on  CON2 
graduation 

3. 

1.-4. 

Dimension- 
less 

INT-1751  estimate 

97 

Unrestricted  CON2 
graduation  rate 

.25 

.2-.4 

Yr"1 

INT-1751  estimate 

98 

CON3  mortality 

.1 

.05-.  15 

Yr"1 

Noble  and  Ronco  1978 

99 

Aspen  restriction  on  CON3 
graduation 

.5 

.3-7 

Dimension- 
less 

INT-1751  estimate 

100 

Conifer  restriction  on 
CON3  graduation 

.3 

.1-.5 

Dimension- 
less 

INT-1751  estimate 
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Param- 
eter FORTRAN 

number         name  Definition  Value  Range  Units  Source 


101 


Unrestricted  CON3  gradua- 
tion rate 


.1-.3 


Yr 


- 1 


INT-1751  estimate 


102 


CON3  graduation 
conversion  to  biomass 


6. 


5.-7. 


Kg/tree 


Long  and  Turner  1975 


103 


104 


105 


106 


107 


108 


109 


110 


111 


112 


113 


114 


Aspen  restriction  on 
conifer  growth 

Aspen  restriction  on 
conifer  growth 

Conifer  restriction  on 
conifer  growth,  mortality  of 
conifer  trees  (#'s)  on  site, 
and  conifer  biomass 
mortality 

Conifer  restriction  on 
conifer  growth  and  unre- 
stricted conifer  growth  rate 

Aspen  restriction  on 
production  of  annuals 

Conifer  restriction  on 
production  of  annuals 

Conifer  restriction  on 
production  of  annuals 

Perennial  restriction  on 
production  of  annuals 

Shrub  and  aspen  sucker 
restrictions  on  production 
of  annuals 

Conifer  seedling  restriction 
on  production  of  annuals 

Conifer  seedling  restriction 
on  production  of  annuals 

Unrestricted  production 
rate  of  annuals 


.5 


3. 


.02 


.07 


2. 


.005 


6. 


.005 


.3-7 


2.-4. 


.01 -.03 


.05-.09 


1.-3. 


.001 -.01 


4.-8. 


.4-.8 


1.-3. 


.001 -.01 


4.-8. 


Dimension-  INT-1751  estimate 
less 

Dimension-  INT-1751  estimate 
less 


Yr 


1 


500 


Yr' 


Dimension- 
less 

Dimension- 
less 

Dimension- 
less 

Dimension- 
less 

Dimension- 
less 


Dimension- 
less 

Dimension- 
less 

kg/ha/yr 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


Bartos  and  Mueggler 
1981 


115 


116 


117 


118 


Aspen  restriction  on 
perennial  production 

Conifer  restriction  on 
perennial  production 

Conifer  restriction  on 
perennial  production 

Annual  restriction  on 
perennial  production 


.007 


6. 


.3-7 


.001-.01 


4.-8. 


.5-.9 


Dimension- 
less 

Dimension- 
less 

Dimension- 
less 

Dimension- 
less 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


119 


Shrub  and  aspen  sucker 
restriction  on  perennial 
production 


.5 


.3-7 


Dimension-  INT-1751  estimate 
less 
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Param- 
eter FORTRAN 
number  name 


Definition 
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Range 


Units 


Source 


120 


121 


122 


123 


124 


125 


126 


127 


128 


129 


130 


131 


132 


133 


134 


135 


136 


137 


138 


Conifer  seedling  restriction 
on  perennial  production 

Conifer  seedling  restriction 
on  perennial  production 

Unrestricted  production 
rate  of  perennials 

Aspen  restriction  on  shrub 
regeneration 

Conifer  restriction  on  shrub 
regeneration 

Shrub  restriction  on  shrub 
regeneration 

Shrub  restriction  on  shrub 
regeneration 

Total  of  aspen  sucker 
restriction  on  shrub 
regeneration 

Total  of  conifer  seedling 
restriction  on  shrub 
regeneration 

Annual  and  perennial 
restriction  on  shrub 
regeneration 

Unrestricted  regeneration 
rate  of  shrub 
regeneration 

Aspen  restriction  on  shrub 
growth 

Conifer  restriction  on  shrub 
growth 

Conifer  restriction  on  shrub 
growth 

Shrub  restriction  on  shrub 
growth  and  shrub  mortality 

Shrub  restriction  on  shrub 
growth  and  unrestricted 
shrub  growth  rate 

Total  of  aspen  sucker 
restriction  on  shrub 
growth 

Total  of  conifer  seedling 
restriction  on  shrub  growth 

Annual  and  perennial 
restriction  on  shrub 
growth 


.007 


.001 -.01 


4.-8. 


3,500 


3.5 


.4 


100 


.001 


.5 


.03 


.01 


2.-6. 


4.-8. 


0-.2 


2.-5. 


.2-.6 


4.-8. 


.6-.9 


75-125 


.4-.6 


0-.01 


.25-1. 


.02-.04 


.15-.25 


.3-7 


0-.1 


.2-.6 


Dimension- 
less 

Dimension- 
less 

kg/ha/yr 


Dimension- 
less 

Dimension- 
less 

Dimension- 
less 

Dimension- 
less 

Dimension- 
less 


Dimension- 
less 


Dimension- 
less 


kg/yr 


Dimension- 
less 

Dimension- 
less 

Dimension- 
less 


Yr 


1 


Yr 


- 1 


INT-1751  estimate 


INT-1751  estimate 


Bartos  and  Mueggler 
1981 

INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


INT-1751  estimate 


Dimension-  INT-1751  estimate 
less 


Dimension-  INT-1751  estimate 
less 

Dimension-  INT-1751  estimate 
less 
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A  deterministic  model  of  succession  in  aspen  forests  was  developed  using 
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allowing  the  parameter  values  to  vary  at  random  within  limits,  was  larger  than  de- 
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in  parameter  values.  These  results  have  indicated  areas  of  needed  research.  The 
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of  eight  regional  experiment  stations  charged  with  providing  scien- 
tific knowledge  to  help  resource  managers  meet  human  needs  and 
protect  forest  and  range  ecosystems. 

The  Intermountain  Station  includes  the  States  of  Montana, 
Idaho,  Utah,  Nevada,  and  western  Wyoming.  About  231  million 
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classified  as  forest  and  rangeland.  These  lands  include  grass- 
lands, deserts,  shrublands,  alpine  areas,  and  well-stocked  forests. 
They  supply  fiber  for  forest  industries;  minerals  for  energy  and  in- 
dustrial development;  and  water  for  domestic  and  industrial  con- 
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Moscow,  Idaho  (in  cooperation  with  the  University  of 
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Reno,  Nevada  (in  cooperation  with  the  University  of 
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